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Abstract. This article gives relations between two types of phase space distributions as- 
sociated to eigenfunctions <pi rj of the Laplacian on a compact hyperbolic surface X-p: 

• Wigner distributions J s . Xr a dWi r . = (Op(a)(f)i rj , 4>i rj }L 2 {x r )i which arise in quantum 
chaos. They are invariant under the wave group. 

• Patterson-Sullivan distributions PSi Tj , which are the residues of the dynamical zeta- 

functions Z(s;a) :— J2 7 \- e -^i I j0 a ( wnere the sum runs over closed geodesies) at 
the poles s = ^ + irj. They are invariant under the geodesic flow. 
We prove that these distributions (when suitably normalized) are asymptotically equal 
as Tj — > oo. We also give exact relations between them. This correspondence gives a new 
relation between classical and quantum dynamics on a hyperbolic surface, and consequently 
a formulation of quantum ergodicity in terms of classical ergodic theory. 
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1. Introduction, statement of results. 




The 


purpose of this article is to relate two kinds of phase space distributions 


which are nat- 



urally attached to the eigenfunctions 0j r of the Laplacian A on a compact hyperbolic surface 
Xp. The first kind are the Wigner distributions Wi Tj G X )/ (5'*Xr)) ()1.1|) of quantum mechan- 
ics. The second kind are what we call Patterson-Sullivan distributions PSi rj G V(S*~Kr)) 
(jl.3|) . We prove that the Patterson-Sullivan distributions have a purely classical mechanical 
definition: they are the residues of classical dynamical zeta functions at poles in the 'critical 
strip'. We also prove that there exists an 'intertwining operator' L r (II. 5|) which transforms 
PS irj — > Wi r . and which induces an asymptotic equality W irj ~ 2 _ ^ 1+2jr ^(— AiTrrj)^ 1 ^ 2 PS irj 
between them. These results indicate a surprisingly close relation between quantum and 
classical dynamics on a compact hyperbolic surface. A similar relation on finite volume 
hyperbolic manifolds of all dimensions is to be expected, but certainly these relations are 
special properties of locally symmetric manifolds. 

To state our results, we introduce some notation. We write G = PSU(1, 1), K = PSO{2) 
and identify the quotient G/K with the hyperbolic disc D. We let F C G denote a co- 
compact discrete group and let X r = T\D denote the associated hyperbolic surface. By 
"phase space" we mean the unit cotangent bundle S**Xr, which may be identified with the 
tangent bundle SX r and also with the quotient T\G. We denote by A = < Ai < A 2 ... the 
spectrum of the laplacian on Xr; with the usual parametrization Aj = Sj(l — Sj) = \ + r| 
(sj = n+^ r i)j we denote by {0j rj }i=o,i,2,... an orthonormal basis of real- valued eigenfunctions: 

A0 irj = -\j<t>i Tj . 

The Wigner distributions (microlocal lifts, microlocal defect measures...) W% r - G P / (5*Xr)) 
are defined by 

(1.1) / a(g)W irj (dg) = (Op(a) ( t) irj , ( j> ir .) LHXr) , aeC°°(S*X r ) 

JS'Xr 
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where Op(a) is a special quantization of a, denned using hyperbolic Fourier analysis (Defini- 
tion E3J). The Wigner distribution Wi rj depends quadratically on (p ir . and has the quantum 
invariance property 

(1.2) (U*Op(a)U t (l) irj ,(f) irj ) = (Op{a)<f) ir ,,<p ir] ), (U t = exp(zVA)); 

hence by Egorov's theorem Wi r . is asymptotically invariant under the action of the geodesic 
flow g l on S"*Xr, in the large energy limit rj — > +00. The Wigner distribution Wi r . is one 
of the principal objects in quantum chaos: it determines the oscillation and concentration of 
the eigenfunction <pi rj in the classical phase space S'*Xr (see $2J). One of the main problems 
in quantum chaos is the quantum unique ergodicity problem of determining which geodesic 
flow invariant probability measures arise as weak* limits of the Wigner distributions (cf. 
|Lin| IRS| IW1 ISh| ISVt IZ2| for a few articles on hyperbolic quotients). 

The Patterson-Sullivan distributions {PSi r } associated to the eigenfunctions {<pi r } (cf. 
Definition 13. 3|) are defined by the expression 

(1.3) PS irj (dg) = PS irj (db',db,dt) := ® \dt\. 

In this definition, T ir is the boundary values of <pi r . in the sense of Helgason (cf . Theorem 
13.11 or |He[ IH].) The parameters (b',b) (b 7^ V) vary in B x B, where B = <9D is the 
boundary of the hyperbolic disc, and t varies in R; (b', b) parametrize the space of oriented 
geodesies, t is the time parameter along geodesies, and the three parameters (&', b, t) are used 
to parametrize the unit tangent bundle ST). These Patterson-Sullivan distributions PSi r . 
are invariant under the geodesic flow g l on SXr, i.e. 

(1.4) (g%PS irj = PS irj . 

The distributions PSi r are also T-invariant (cf. Proposition l3.3|) . Hence they define geodesic- 
flow invariant distributions on 5*Xr. Phase space distributions of this kind were associated to 
ground state eigenfunctions of certain infinite area hyperbolic surfaces by by S. J. Patterson 
|FatUl iPatT] . and were studied further by D. Sullivan |Sull ISu2j (see also jNj). Ground 
state Patterson-Sullivan distributions are positive measures, but our analogues for higher 
eigenfunctions on compact (or finite area) hyperbolic surfaces are not measures. To our 
knowledge, they have not been studied for higher eigenfunctions before. 

Both families (Wi r .) and (PSi rj ) are T- invariant bilinear forms in the eigenfunctions <fii r . 
with values in distributions on SXp. But they possess different invariance properties: the 
former are invariant under the quantum dynamics (the wave group) while the latter are 
invariant by the classical evolution (the geodesic flow). The motivating problem in this 
article is to determine how they are related. 

Our first result is that, after suitable normalization of PSi rj , the Wigner and the Patterson- 
Sullivan distributions are asymptotically the same in the semi-classical limit rj — > +00. The 
exact relation between involves the operator L r : C^(G) — > C°°(G) defined by 

(1.5) L r a{g)= [ (1 + u 2 )-^ +ir) a{gn u )du 

Jr 
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which, we will see, mediates between the classical and quantum pictures. Here, n u = 
1 u \ 

acts on the right as the horocycle flow. We further introduce a cutoff function 

ij 

X G Co°(D) which is a smooth replacement for the characteristic function of a fundamental 
domain for F (called a 'smooth fundamental domain cutoff', see Definition IH.2j) . 

Theorem 1.1. For any a e C°°(r\G) we have the exact formula 

(Op(a)0 4rj ,0„,) = 2 -( 1+2 ^) / (L rjX a)dPS irj , 
and the asymptotic formula 

_ a(g)W ir .(dg) = 2^ 1+2 ^)(-4mr J )- 1/2 f a(g)PS irj (dg) + O^). 

It follows that the Wigner distributions are equivalent to the Patterson- Sullivan distribu- 
tions in the study of quantum ergodicity. The operators L r in a sense intertwine classical 
and quantum dynamics. We note that, although the Wigner distributions were defined by 
using the special hyperbolic pseudodifferential calculus Op, any other choice of Op will pro- 
duce asymptotically equivalent Wigner distributions and hence Theorem II .11 is stable under 
change of quantization. 

When a is an automorphic eigenfunction, i.e. a joint eigenfunction of the Casimir operator 
Q and the generator W of K, we can evaluate the first expression in Theorem 11.11 to obtain 
a very concrete relation: 

Theorem 1.2. If a is an eigenfunction of Casimir parameter r and weight m in the con- 
tinuous series, we have: 
(i) 

{Op{a)<j> ir .^ irj ) = ^(l+irj) (j nG o-dPS irj ) +/#(i + ir i ) (j r \ G X + cr dPS tr . 

where \i c mr and \i c ^ are defined in (|5.b|) : X + denotes the vector field generating the horocycle 
flow. 

(ii) If a = ip m is a lowest weight vector in the holomorphic discrete series, we have 
(Op(ip m )(j) ir] ,(i) ir] ) = n m (^ + irj) (^J ^ adPSi r 

where \x m is defined in \5. 7| ). 

These exact formulae are based on the identity (cf. Proposition 16. 4|) . 

(1.6) / (L r . X <r)dPS irj = [ (l + u 2 )- {1 i +irj) Ips ir .(o-)(u)du, 

Jsr> Jn 3 

where Ips iT] '■ C°°(T\G) — > C(R) is the operator defined by 

(1-7) I PSl (o-)(u):= I a(gn u )PS irj (dg). 

Jt\g 
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When a is a joint eigenfunction of the Casimir operator f2 and of the generator W of the 
maximal compact subgroup K, the function Ips ir (a)(u) is a special function of hyperge- 
ometric type depending on r and the eigenvalue parameters of a (cf. ^21 for a review of 
the representation theory of L 2 (T\G)). The integral on the right side of (jl.6|) can then be 
evaluated to give the explicit formulae of Theorem 11.21 

Theorems 11.11 and 11.21 admit generalizations to off-diagonal matrix elements, which relate 
off-diagonal Wigner and Patterson-Sullivan distributions by operators L r>r /. These gener- 
alizations will be explored in our subsequent article |AZj . It is obvious that there cannot 
exist any standard kind of intertwining between the quantum dynamics (which has a discrete 
spectrum) and the classical dynamics (which has a continuous spectrum), but the operators 
L r y establish a kind of intertwining on the level of invariant distributions. 

Our next result gives a purely classical dynamical interpretation of the Patterson-Sullivan 
distributions in terms of closed geodesies. Given a e C°°(SXt), we define two closely related 
dynamical zeta-functions 

[ (<) Z(s; a) := £ 7 a) , s > 1), 

(1.8) 

where the sum runs over all closed orbits, and 70 is the primitive closed orbit traced out by 
7. The sum converges absolutely for die s > 1. 

Theorem 1.3. Let a be a real analytic function on the unit tangent bundle. Then Z(s; a) and 
^2(3; a>) admit a meromorphic extension to C. The poles in the critical strip < 3te s < 1, 
appear at s = 1/2 + ir, where as above 1/4 + r 2 is an eigenvalue of A. 
The residue is 

— -.2 2 

].r-—r 

where {dPSi Tj } are the Patters on- Sullivan distributions associated to an orthonormal eigen- 
basis {4>irj}, and /xo is defined in ()5.6|) . 

The normalisation factor Ho(\ + irj) is such that /-tod + ^j) / UPSi rj = 1, for a Patterson- 
Sullivan distribution associated to an eigenfunction normalized in L 2 . Note that it does not 
depend on the group T. 

We give two quite different proofs of this theorem, each one providing a different perspec- 
tive on the correspondence between Wi Tj and PSi rj . The first proof uses the thermodynamic 
formalism. It is known that this formalism can be used to locate the zeros of Selberg's zeta 
function |Polj . In ^3 we apply similar ideas to prove that Z2 admits a meromorphic contin- 
uation to 3?e s > 0, and we describe its poles and residues in terms of "Ruelle resonances". 
In particular, Patterson-Sullivan distributions arise as the residues. We use the methods 
developed by Rugh in Rugh92 1 for real-analytic situations. The techniques are based on the 
Anosov property of the geodesic flow, and apply in variable curvature. However, the relation 
between Wigner and Patterson-Sullivan distributions is special to constant curvature. 

The second proof, given in $HJ uses representation theory and the generalized Selberg trace 
formula of [Z] to prove the meromorphic continuation of Z to all of C and to determine the 
poles and residues. Like Theorems 11.11 - 11.21 the trace formula establishes an exact relation 
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between the Wigner distributions (which appear on the 'spectral side' of the trace formula) 
and the geodesic periods fa on the 'sum over T' side. No such formula can be expected in 
variable curvature, and the methods are specific to hyperbolic surfaces. 

In conclusion, the results of this paper develop to a new level the strangely close relation 
between classical and quantum dynamics on hyperbolic surfaces. On the level of eigenvalues 
and lengths of closed geodesies, this close relation is evident from the Selberg trace formula 
(cf. 3EJ). As is well-known, the Selberg trace formula on a compact hyperbolic manifold is 
a special case of the general wave trace formula on a compact Riemannian manifold where 
the leading order approximation is exact. The exactness of this stationary phase formula 
is somewhat analogous to the exact stationary phase formula of Duistermaat-Heckman for 
certain oscillatory integrals, but to our knowledge no rigorous link between these exact 
formulae is known. An alternative explanation of the close relation between classical and 
quantum dynamics was suggested by V. Guillemin in [U|, who made a formal application of 
the Lefschetz formula to the action of the geodesic flow on a non-elliptic complex. The trace 
on chains gave the logarithmic derivative of the (Ruelle) zeta function, while the trace on 
homology gave the spectral side of the Selberg trace formula. For later developments in this 
direction (by C. Denninger, A. Deitmar, U. Bunke, M. Olbrich and others) we refer to jj]. 

This paper develops the close relation on the level of eigenfunctions and invariant distri- 
butions rather than just eigenvalues and lengths of closed geodesies. In future work |AZj . 
the intertwining of classical and quantum mechanics will be developed further so that it 
applies not just to traces (diagonal matrix elements) but to the full dynamics. It is hoped 
that it will have applications related to the question of quantum unique ergodicity. It would 
also be interesting to relate our constructions to generalizations of the non-elliptic Lefschetz 
formulae and to other representation theoretic ones in |SVt IWj. 

Acknowledgements This work was begun while the first author was visiting Johns Hopkins 
University as part of the NSF focussed research grant # FRG 0354386. Much of it was 
written at the Time at Work program of the Institut Henri Poincare in June, 2005. 



2. Background 

Hyperbolic surfaces are uniformized by the hyperbolic plane H or disc D. In the disc 
model D = {z e C, \z\ < 1}, the hyperbolic metric has the form 

ds 2 - 



1- z 



2\2 ' 



The group of orientation-preserving isometries can be identified with PSU(1, 1) acting by 
Moebius transformations; the stabilizer of is K ~ 5*0(2) and thus we will often identify D 
with SU(1, 1)/K. Computations are sometimes simpler in the H model, where the isometry 
group is PSX(2, R). We therefore use the general notation G for the isometry group, and 
G/K for the hyperbolic plane, leaving it to the reader and the context to decide whether 
G = PSU(1, 1) or G = PSL(2, R). 

In hyperbolic polar coordinates centered at the origin 0, the Laplacian is the operator 

d 2 8 Id 2 

A = — — + cothr— — h 



dr 2 dr sinh 2 r d6 2 
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The distance on D induced by the Riemannian metric will be denoted g?d- We denote the 
volume form by dVol(z). 

Let r C G be a co-compact discrete subgroup, and let us consider the automorphic 
eigenvalue problem on G/K: 

( A0=-A0, 

(2-1) 

I 0(72) = 4>{z) for all 7 G Y and for all z. 

In other words, we study the eigenfunctions of the Laplacian on the compact surface 
Xr — Y\G/K. Following standard notation (e.g. [310]), the eigenvalue can be written in 
the form A = A r = 7 + r 2 and also A = A s = s(l — s) where s = | + ir. 

Notational remarks 

(i) We denote by {Xj = -g + r 2 } the set of eigenvalues repeated according to multiplic- 
ity, and (in a somewhat abusive manner) we denote a corresponding orthonormal basis of 
eigenfunctions by {4>i rj }- 

(ii) We follow the notational conventions used in jN] and [0] . which differ from those used 
in jHj, where the metric is multiplied by a factor 4. We caution that [LIE] use the latter 
conventions, and there the parameter s is defined so that 4A = (s — l)(s + 1) and so that 
s = 2ir. 

2.1. Unit tangent bundle and space of geodesies. We denote by B = {z G C, \z\ = 1} 

the boundary at infinity of D. The unit tangent bundle SD of the hyperbolic disc D is 
by definition the manifold of unit vectors in the tangent bundle TD with respect to the 
hyperbolic metric. We may, and will, identify SD with the unit cosphere bundle S*D by 
means of the metric. We will make a number of further identifications: 

• SD = PSU(1,1). This comes from the fact that PSU(1,1) acts freely and tran- 
sitively on SD. We identify a unit tangent vector (z, v) with a group element g if 
g- (i, (0, 1)) = (z, v). Similarly, if we work with the upper half plane model H, we have 
SH = PSL(2,R). We identify SD, SH, PSU(1, 1), and PSL(2,R). In general, we 
work with the model which simplifies the calculations best. According to a previous 
remark, SD, PSU(1, 1) and PSL(2, R) will often be designated by the letter G. 

• STJ = D x B. Here, we identify (z, b) G D x B with the unit tangent vector (z,v), 
where v G S Z D is the vector tangent to the unique geodesic through z ending at b. 

The geodesic flow g t on SD is defined by g t (z,v) = j' v (t)) where j v (t) is the unit 

speed geodesic with initial value (z,v). The space of geodesies is the quotient of SD by the 
action of g*. Each geodesic has a forward endpoint b and a backward endpoint b' in B, hence 
the space of geodesies of D may be identified with B x B \ A, where A denotes the diagonal 
in B x B: To (b', b) G B x B \ A there corresponds a unique geodesic 7^ whose forward 
endpoint at infinity equals b and whose backward endpoint equals b' . 

We then have the identification 

SD = (B x B \ A) x R. 

The choice of time parameter is defined - for instance - as follows: The point (6', 6, 0) is by 
definition the closest point to on 7^ and (b',b,t) denotes the point t units from (b, b',0) 
in signed distance towards b. 



8 NALINI ANANTHARAMAN AND STEVE ZELDITCH 

2.2. Non-Euclidean Fourier analysis. Following [H], we denote by (z, b) the signed 
distance to of the horocycle through the points z G D, b G B. Equivalently, 

e{z,b) = y=w = p » {z > b) > 

where -Pd(^, b) is the Poisson kernel of the unit disc. (We caution again that e^'^ is written 
e 2 ( z > b ) j n [h | IZ2j ). We denote Lebesgue measure on B by \db\, so that the harmonic measure 
issued from is given by Pn(z, b)\db\. A basic identity (cf. [H]) is that 

(2.2) (g-z,g-b) = (z,b) + (g-0,g-b), 
which implies 

(2.3) P D {gz,gb)\d(gb)\ = P D (z,b)\db\. 

The functions e^ +ir )( z > b ) are hyperbolic analogues of Euclidean plane waves e*^'^ and are 
called non-Euclidean plane waves in jH] . The non-Euclidean Fourier transform is defined by 

Fu(r,b) = [ e^- ir){z ' b) u(z)dVol(z). 

The hyperbolic Fourier inversion formula is given by 

u{z) = [ [ e { 5 +ir )^^u(r,6)rtanh(27rr)rfr|d6|. 



B JR. 



As in |Z3j . we define the hyperbolic calculus of pseudo-differential operators Op(a) on D by 

Op(a)e^ +ir){z > b) = a(z,b,r)e^ +ir ^ z ' b) . 

We assume that the complete symbol a is a polyhomogeneous function of r in the classical 
sense that 

oo 
3=0 

for some m (called its order). By asymptotics is meant that 

R 

a{z, b,r)-J2 °i(^> b)r~ j+m G S^^ 1 

3=0 

where a G S k if sup(l + ry~ k \D°D^Dj,(z)\ < +oo for all compact set and for all a,j3,j. 

The non-Euclidean Fourier inversion formula then extends the definition of Op(a) to 
C?°(D): 



Op(a)u(z) = / / a(z,b,r)e^ 2 +ir ^Fu(r, 6)rtanh(27rr)dr|d6|. 



B JR. 



A key property of Op is that Op(a) commutes with the action of an element 7 G G 
(Tryu(z) = u{^z)) if and only if 0(7,2,76, r) = a(z, b, r). T-equivariant pseudodifferential 
operators then define operators on the quotient Xp. This will be seen more clearly when we 
discuss Helgason's representation formula for eigenfunctions. 
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2.3. Dynamics and group theory of G = SL(2, R). We recall the group theoretic point 
of view towards the geodesic and horocycle flows on SXp. As above, it is equivalent to work 
with G = PSU(1, 1) or G = SX(2,R); we choose the latter. Our notation follows [El [Zj, 
save for the normalization of the metric. The generators of sl(2, R) are denoted by 



H=\ , V=\ , W 






We denote the associated one parameter subgroups by A,A_,K. We denote the rais- 
ing/lowering operators for i^-weights by 

(2.4) E + = H + iV, E~ = H- iV. 

The Casimir operator is then given by 4 Q = H 2 + V 2 — W 2 ; on K-invariant functions, the 
Casimir operator acts as the laplacian A. We also put 









( ° 




: 








: 













and denote the associated subgroups by N, N_. 

In the identification SD = PSL(2, R) the geodesic flow is given by the right action of the 
group of diagonal matrices, A: g t (g) = ga t where 




at 



a 

By a slight abuse of notation, we sometimes write a for ( | . The action of the 

a' 1 

geodesic flow is closely related to that of the horocycle flow (h u ) ue n, defined by the right 
action of N, in other words by h u (g) = gn u where 




n u = 

The closed orbits of the geodesic flow g l on T\G are denoted {7} and are in one-to-one 
correspondence with the conjugacy classes of hyperbolic elements of T. We denote by G 7 , 
resp. r 7 , the centralizer of 7 in G, resp. T. T 7 is generated by an element 70 which is called 
a primitive hyperbolic geodesic. The length of 7 is denoted L 1 > and means that 7 is 
conjugate, in G, to 



(2.5) a 



7 



e W2 
y e"W 2 

If 7 = 7q where 70 is primitive, then L 7o is the primitive length of the closed geodesic 7. 
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2.4. Representation theory of G and spectral theory of A. Let us recall some basic 
facts about the representation theory of L 2 (r\G) in the case where the quotient is compact 

(cf. mu)- 

In the compact case, we have the decomposition into irreducibles, 

5 oo oo oo 

L 2 (r\G) = 0C ! ,©0P irj © Mm)D+© A*r(m)P-, 

j=l j=0 m=2, m even m=2,m even 

where C Sj denotes the complementary series representation, resp. V Sj denotes the unitary 
principal series representation, in which — Q equals Sj(l — Sj) = f +?"|- hi the complementary 
series case, irj G R while in the principal series case irj G The irreducibles are indexed 
by their i^-invariant vectors .}, which is assumed to be the given orthonormal basis 
of A-eigenfunctions. Thus, the multiplicity of Vi rj is the same as the multiplicity of the 
corresponding eigenvalue of A. 

Further, T>^ denotes the holomorphic (resp. anti-holomorphic) discrete series representa- 
tion with lowest (resp. highest) weight m, and fir( m ) denotes its multiplicity; it depends 
only on the genus of X r . We denote by ip m j (j = 1, . . . , fj,r(m)) a choice of orthonormal basis 
of the lowest weight vectors of /zr(m)D^ and write jjr( m )^m = ®t=i j accordingly. 

We will also use the notations Ci rj} Vi rj and T> m j for the orthogonal projection operators 
of L 2 (r\G) onto these subspaces. Thus, for / G L 2 we write 

(2.6) / = »,,.(./•.) + £>*.,(/) + £ ■n;„Jf). 

By an automorphic (r, m)-eigenfunction, we mean a T-invariant joint eigenfunction 



(2.7) 



Wa Tm = ima T _ m . 



of the Casimir Q and the generator W of K = 50(2). 

We recall that the principal series Vf r are realized on the Hilbert space L 2 (R) by the 
action 

^/w = l-^l— /(^). 

The unique normalized i^-invariant vector of Vi Tj is a constant multiple of 

U, (x) = (l + x 2 )-^ +ir \ 
The complementary series representations are realized on L 2 (R, B) with inner product 



and with action 




W.f> = J 

Jrxr \ x y\ 



f(x) = \-bx + d\- l - 2u f(- aX ° 



—bx + d' 
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When asymptotics as \rj\ — > oo are involved, we may ignore the complementary series 
representations and therefore do not discuss them in detail. 

Let C + = {z E C : Qz > 0}. We recall (see jK], §2.6) that X>+ is realized on the Hilbert 
space 

= {/ holomorphic on C+, / \f(z)\ 2 y m ~ 2 dxdy < oo} 

Jc+ 



with the action 




v + m \ \f(z) = (-bz + dy™f(^L). 

The lowest weight vector of T*^ in this realization is is (z + i)~ m . 

We note that the .fT-weights in all irreducibles are even. Lowest weight vectors of 
correspond to (holomorphic) automorphic forms of weight m for V in the classical sense of 
holomorphic functions on C + satisfying 



f(i-z) = (cz + d) m f(z), 7 = , 7 er. 




Thus, / is invariant under the action of T. A holomorphic form of weight m defines a 
holomorphic differential of type f(z)(dz)^ (cf. |Sa2j ). Forms of weight n in Vi r ,C u ,V^. 
always correspond to differentials of type (dz)*. Forms of odd weight do not occur in 
L 2 {T\PSL{2,R)). 

We now consider the action of A, i.e. the geodesic flow, in each irreducible. 

Proposition 2.1. The right action of A, i.e. the geodesic flow g l , has two invariant sub- 
spaces in each irreducible Ci r ,Vi r , namely the cyclic subspace generated by the weight zero 
vector (fi ir , and that generated by X + (p ir . The action of A is irreducible in Z>*. 

Proof. In the principal series we have 

/ e'/ 2 \ 
Vl fix) = S^f&x). 

\ e-*/ 2 / 

The subspaces L 2 (R + ), L 2 (R_) are invariant. The generalized eigenfunctions (eigendistri- 
butions) . The action is irreducible in each subspace. 

In the discrete series we have 

e 1 ' 2 

Z£ I I /(*) = e m4/2 /(e^) 

e-*/ 2 

The lowest weight vector is cyclic for the action of A. 



□ 
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2.5. Time reversibility. Time reversal refers to the involution on the unit cosphere bundle 
defined by i(x,£) = (x, — £). Under the identification T\G ~ S*Xr, the time reversal map 

/0 1\ 

takes the form Tg — > Twgw where w = I is the Weyl element. For a G A one has 

V- 1 °/ 

waw = a 1 . 

We say that a distribution is time-reversible if l*T = T. The distributions of concern in 
this article all have the property of time-reversibility, originating in the fact that A is a real 
operator and hence commutes with complex conjugation. This motivates the decomposition 
of Vir = Vf r © V~ into 'even' and 'odd' subspaces. 

Proposition 2.2. We have: 

• Each principal (or complementary) series irreducible contains a one- dimensional 
space of A-invariant and time-reversal invariant distributions. In the realization on 
L 2 (R) 7 it is spanned by £ r (x) = \x\~ ( -2 +lT '>\ 

• There exists a unique (up to scalars) A-invariant time-reversal invariant distribution 
in when m = 0(modA) and there exists no time reversal invariant distribution 
when m = 2{modA). In the realization on 7i+ ; it is z~' m l 2 . Similarly for X>~. 



Proof, (i) The complementary and principal series 

Each principal (or complementary) series irreducible contains a two-dimensional space 

of A-invariant distributions. In the model on £ 2 (R) a basis is given by 
Indeed, A invariance is equivalent to 

e-« 1 > +i % r (e t x)=Z ir (x). 

Setting x = ±1 we find that 

are invariant distributions supported on R±. 
The time reversal operator is given by 

I ° l \ -1 

(2.9) V ir f(x) = \ x \-^f(— ). 

\ -1 / x 
Hence, time reversal invariance is equivalent to 

/(— ) = \x\ 1+2ir f(x). 

Ob 

Under time reversal 

Vi r (w)X + 2 = \X\ X_ = X_ 2 

Hence the unique time reversal invariant distribution is 

fc r = \x\~^ +ir \ 



(ii) The discrete series 
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Each holomorphic (or anti-holomorphic) discrete series irreducible contains a unique 
(up to scalar multiple) A- invariant distribution z~ m l 2 . Indeed, to solve 

/ e'/ 2 \ 

V e-*/ 2 / 
we put z = e %e and obtain 

and the only holomorphic solution is z~ m l 2 . 

In the holomorphic discrete series, the time reversal operator is given by 

We observe that 

z -m/2 ig 

time-reversal invariant when m = 0(mod4) and is anti-invariant 

when m = 2mod4. 

The anti-holomorphic discrete series is similar (by taking complex conjugates). 

□ 

Definition 2.1. We denote the time reversal and geodesic flow invariant distribution in 
V(T\G) n Vir j: normalized so that (0^,2^) = 1, by E irj . We denote by ■ the 
time reversal and geodesic flow invariant distribution in D'(T\G) fl T>^, normalized so that 
{ip m j,E m j) = 1, where ||Vwll — 1- Here, we assume m = 0(4). 

The following is the main application of the representation theory. By the above normal- 
ization, all denominators equal one, but we leave them in to emphasize the normalization. 

Proposition 2.3. Let dv denote a time-reversal invariant and geodesic flow invariant dis- 
tribution on T\G. Let f eC°°(T\G). Then: 

I V^ 00 ( P m,j/' 5 m,j) / / J \ 

+ ^±^2,^0(4)^=1 ^~^~) KWm^aV). 

Proof. Since <pi r is cyclic in Vi r , any element / in this space may be expressed in the form 
J R /(t)0i r o g t dt. If we pair with the invariant distribution dv we obtain J R f(t)dt ((pi r ,du). 
On the other hand, if we pair / with E ir we obtain f R f(t)dt ((f) ir ,E ir ) . Similarly in the 
discrete series. The statement follows immediately. 

□ 

To apply the Proposition, we need to understand convergence of the series and hence 
to have bounds on (Vi r (f), Hj r ) and (X>^ •/, •) when the denominator is normalized to 
equal one. Since the complementary series sum is finite, it is not necessary to analyze these 
terms. The following proposition shows that the distributions are of order one. Here, we say 
that a distribution T has order r if (T,f) < \\f\\w s where W S (T\G) is the Sobolev space 
of functions with s derivatives in L? . The proposition also controls the dependence of the 
norms in the Casimir parameters ir, m. 
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Proposition 2.4. We have: 

• \(PUf),E r )\ <Cr-V 2 \\V ir (f)\\ w i; 
. \(ViJ,~ m ,)\<Cm-^\\ViJ\\ m; 

Proof. We prove the results by conjugating to the models above. 
We begin with the continuous series and let 

U ir : L 2 (T\G) -> L 2 (R) 

be the unitary intertwining operator from Vi r C L 2 (T\G) to its realization in L 2 (R). Thus, 
Uir'B.ir = & r up to the normalizing constant. 

To determine the normalizing constant, we recall (see [Zj, p. 59) that 

(UiAirMir^ir) = / R (l + X 2 )^ + ^ \ X \~^dx 



2f °°(l + x 



= 2 J oo (a;- 1 + x)-(^)f 
= 2B(j(J + ir),|(j + <r)):=2^^ 
Here, B(x,y) = r?^x ^ s ^ ne beta-function. From the asymptotics (cf. |GR| 8.328) 

(2.10) F{x + iy) ~ V27re"5 |tfl |3/| a! "3 (|y| -> oo) 

of the T-function along vertical lines in C, it follows that 

TV I -L Ie^ 

(A-)" 1 := r ; 4 ^ 2 ; ~ Cr-^, (r _> oo). 

Next we consider the order of £j r as a distribution in the model. We may break up each 
function in L 2 (R) into its even and odd parts with respect to time reversal invariance, and 
then we only need to consider (/,£«•) for a time reversal invariant /. Let x+i x ) Co° with 
X+ = 1 for |x| < | and for |x| > 2 and with the property that x+( x ) + X+( ) — 1- Then 
ifi&r) — ((x+ + X+( Cir) an d (denoting the time reversal (|2.9|) operator by T) 

(x+(=r)/,&-> = (Tx+{=±)f,n ir ) 

= {{Tx+i^T-^fMr) 

= (X+fitir)- 

Hence we only need to estimate the x+ integral. We write x~ l l 2+tT = _ x J 2 +ir ~Ex xl ^ 2+%r ano - 
integrate by parts. The result is bounded by C(l + r)~ 1 (||/|| i 2 + \\Vi r (X_)f\\ L 2). Here, we 
use that X + is represented by j-. 



PATTERSON-SULLIVAN DISTRIBUTIONS AND QUANTUM ERGODICITY 15 

It follows that for / e C°°(T\G), 

\$^\ = \(V ir (f),E ir )\ = MUrVUf),^ 

— Pr\\-^L(ir'Pir(f)\\L 2 (R) 

= C^(l + r)- 1 ||X_P ir (/)|U 2( r\ G ) < r-Va \\V ir (f)\\ w i { r\ G) . 
We now consider the discrete series. The normalizing constant is now 

IK z ~r^ II JC+ 

To calculate the constant, we use the isometry 



T m : H+ -> 2 (D,^ m ), T m /(«,) = / f-<^4) (— ^ 

\ w — 1 J \w — 1 



where (9 2 (D, cZv m ) are the holomorphic functions on the unit disc which are L 2 with respect 

. 4 f-i \ \2\m_dwdw_ 

/ m — im \ L \ UJ \ ) 



to the measure dv m = ^(1 - |w| 2 ) m totf ( cf. |Lj IX §3). 



We have T m vb m = 1. Note that 1 is not normalized to have L 2 norm equal to one. It 
follows that 

/-+ / \ 4 f ( .w + l\~ m/2 ( -2i \ m 2 dwdw 



4 m ||(z + z)- m || y D V w-iy v^- 1 / (l — |tw 



2\2 



We write w = re and observe that the angular integral equals an r-independent constant 
times 



l+ re «V m/ V 2z \ m ,„ f fl + rzY m/2 ( -2i \ m dz 

d6 = / — = 2vri -2z 



i \ 1 — re™ / \ 1 — re™ J Ji z i =1 \l — rz J \ 1 — rz J z 

since (yz^r) m 2 {f^zi)™ 1S holomorphic in \z\ < 1 for r < 1. It follows that 



since the L 2 -norm of T m ^ m = 1 equals £ (/ D (l - M 2 ) m (iq^w) and /d( X 
equals —^—r. 

^ m— 1 

We then need to estimate 

(-+,;>= / /(z)z- m/2 y m - 2 ^. 



1/2 

'2\m dwdw 



\W\ 



(1-|H 2 ) 2 



As above, we let x+ be a radial function with compact support in R + and with x+( z ) + 
X+(—) = 1- By unitary of time reversal, we again have 

and thus it suffices to estimate the x+ integral. We note that for m > 2, z~~ m l 2 = 
_j_d_ z - m /2 + i and that z -m/2+i e L 2 ( | z | < X ,y md -^L). The operator ± = is 
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skew symmetric with respect to the inner product. Partial integration gives the bound 
i-to/2 H/'IU 2 ) nence after normalizing we have 



<Cm- 1 / 2 ||/|| £a + |K(X + )/|| £a ). 

□ 



Remark 2.1. The paper |A-Pj studies related estimates in the discrete series from a different 
point of view. 

A nice simplification is that the series {X + <pi rk } automatically has zero integrals against 
a time reversal invariant distribution: 

Lemma 2.5. IfT G T>'(T\G) is time-reversible, then (X + <fti rk ,T) = for all k. 
Proof. We have 

(X +( j) ir ,T) = (X +( j) ir ,L*T) 
= (i*{X +( f> ir ),T) 
= -(X + <f> ir ,T). 

□ 

3. Patterson-Sullivan distributions and microlocal lifts 

3.1. Patterson-Sullivan distributions. Let us first recall Helgason's fundamental result 
about eigenf unctions of the Laplacian on D. In the following theorem, <j) is any solution of 
Acf) = —\(f) (A = \ + r 2 where A, r 6 C). The function (ft, defined on D, is not necessarily 
automorphic. One says that has exponential growth if there exists C > such that 
\(p(z)\ < Ce Cd »(°^ for all z. 

Theorem 3.1. (^Hj 7 Th eorems 4-3 and 4-29; see also |Hej ) Let <fi be an eigenfunction with 
exponential growth, for the eigenvalue A = \+r 2 G C. Then there exists a unique distribution 
G V'(B) such that 



f e ^ +ir){z ^T^db) } 
Jb 



for all z G D 



The theorem extends the classical representation theorem for bounded harmonic functions 
to the case of arbitrary eigenvalues. Note that the kernel e^ +ir ^ z ' b ' ) that appears in the 
representation theorem for eigenf unctions for the eigenvalue A r is the generalized Poisson 

kernel, P D 2 (z,b). The distribution is called the (radial) boundary value of and is 
defined as follows: expand the eigenfunction into the "Fourier series" , 

(3.1) <K^) = $> n $ rin (z), 

nGZ 

in terms of the generalized spherical functions $ r n defined by 

(3.2) S +ir){z ' b) = ®r,n( z )b n ' b G B. 

n<=Z 
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The boundary value of is the distribution on B defined by 
(3.3) T^db) = J2 a nb n \db\. 

ngZ 

As observed in |Z2| , when 0j rj is a T- invariant eigenf unction, the boundary values Ti rj (db) 
have the following invariance property: 

fan M = K s (*) => e^ +ir ^ z ^T irj (d 7 6) = e^ r ^T ir . (db) 

(3.4) 

=}► T irj .(d 7 6) = e -(2+^)<ro,r6> T . r .( d6 ) 

This follows from the uniqueness of the Helgason representation (|3.1j) and by the identities 
d23)-(Q. Hence there exist e irj G V'(T\PSL(2, R)) such that 

(3.5) e^^T* . = e*,.. 6)P(2, b)db. 

The distribution ej rj is horocyclic- invariant and T- invariant. It may be expanded in a K- 
Fourier series, 

neZ 

and it is easily seen (cf. |Z2j ) that 0i rj -,o = and that ^i r . iTl is obtained by applying the 
nth normalized (i.e. unitary) raising/ lowering to <pi rj . More precisely, one applies (E ± ) n 
(|2.4j) and multiplies by the normalizing factor fyir^n — (2ir-+i±2n)-(2ir-+i±2) - ^ ne regularity 
of these distributions was recently studied in |FF| ICoj . 

If we set z = 0, the i^-Fourier series and 5-Fourier series coincide and we get 

(3.6) T ir ,{db) = e trj (0,b)db = J2^,±n ((^ ± ) n 0,, J (O)) b n db. 

neZ 

Rather than estimating the regularity of T irj (db) using ()3.6|) . which would take us too far 
afield, we will quote some estimates of Otal [0] which suffice ( and indeed are better than 
necessary) for our applications. Roughly, the say that T ir .(db) is the derivative of a Holder 
continuous function Fi rj . Since its zeroth Fourier coefficient is non-zero, Ti r .{db) is not 
literally the derivative of a periodic function, but it is the derivative of a function F ir . on R 
satisfying Fi Tj [9 + 2ir) = Fi T] iff) + Cj for all 6R. We follow Otal in calling such a function 
27r-periodic. 

For < 5 < 1 we say that a 2vr-periodic function F : R -> C is 8- Holder if \F(9)-F(9 r ) | < 
C\6 — 9'\ s . The smallest constant is denoted \\F\\s and denotes the Banach space of 5- 
Holder functions, up to additive constants. 

Theorem 3.2. (]Uj Proposition 4) Suppose that a = 5 + it and that <fi is an eigenfunction 
of eigenvalue a{\ — a) satisfying ||V0||oo < oo. Then its Helgason boundary value is the 
derivative of a 5-Hdlder function. 

This estimate holds for arbitrary bounded eigenfunctions. In our case, it says that Ti r . is 
the derivative of a Holder function, of Holder exponent 6 = \ if Xj > j. Otal's proof also 
shows that the Holder norm is bounded by a power of rj. Related results can be found in 

pmin miFFi inoi. 
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We now introduce a "Patterson- Sullivan" distribution associated to each automorphic 
eigenfunction. We denote Ao = < Ai < ... the spectrum of the laplacian on Xr (Xj = 
(4>i r ) is a given orthonormal basis of eigenf unctions, and their boundary values are denoted 

(<•)• 



Remark 3.1. We assume that these eigenfunctions are real to obtain time reversal invariant 
distributions. Aside from that, our results are valid for complex eigenfunctions with slight 
modifications. 

Definition 3.1. The Patterson-Sullivan distribution associated to a real eigenfunction 0j r . 
is the distribution on B x B \ A defined by 

T ir .{db)T ir .{db') 
ps irj {db , db) := \ h _ hl ^ 2irj 

If <pi r - is automorphic, it is easy to check that pSi T . is invariant under the diagonal action 
ofT: 

Proposition 3.3. Suppose that cj)i Tj is T-invariant, and let T irj denote its radial boundary 
values. Then the distribution on B x B \ A defined by 

. . , T ir .(db)T ir \db') 
ps irj (db , db) := ^_ y|1+ J 2trj 

is T-invariant and time reversal invariant. 
Proof. It follows from ()3.4|) that 

(3.7) T irj (d 1 b)T ir .(d>yb') = e-^ + ^)^- ^ f> >e-(H^)(7-o,Tfe'> r ^ (rf6)T ^ ( ^). 

We will also need the following identities (cf. |H] (1.3.2)): 

|7(x) -i{y)\ = \j'(x)\3\~/'(y)\?\x-y\ 

(3.8) 

1 - | 7 ( x )|2 = | 7 '(a;)|(l- \x\ 2 ). 
for every x, y G D U B, 7 G T. Hence for b G B and 7 G T, we have 

(3.9) | 7(0 )_ 7 ( & )|2 = |y( 6 )|(i_| 7 (o)| 2 ). 
Furthermore, 

(3.10) | 7& _ 76 '|2 = e -[<rO,T6>+<7-o 1 r6'>]| 6 _ 6 '|2_ 

Raising (|3.1()jl to the power | + irj, taking the ratio with ()3.7|) and simplifying completes 
the proof of r-invariance. 

Time-reversal invariance is invariance under b <^=^ b', which is obvious from the formula. 

□ 

We now construct from the distribution pSi Tj a geodesic flow invariant distribution on 
ST) as follows. As reviewed in £0 the unit tangent bundle can be identified with 
(BxB\A)xR: the set B x B \ A represents the set of oriented geodesies, and R gives 
the time parameter along geodesies. We then define the Radon transform: 

(3.11) K:C (SD)^C (BxB\A), by Kf(b', b) = [ fdt. 
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Further, we need to define special cutoffs which have the property that 

(3.12) / fdVol(z) = [ xfdVol(z) 

Jv Jd 

for any / G C(r\D), where T> is a fundamental domain for T in T>. In other words, \ is a 
smooth replacement for the characteristic function of T>. 

Definition 3.2. We say that x £ Co°(D) is a smooth fundamental domain cutoff if 

= 1. 

We then make the basic definitions: 

Definition 3.3. (1) On ST) we define the Patterson-Sullivan distribution PS ir . eV'(ST)) 
by: 

PS irj (db', db, dt) = psir 3 (db' , db)\dt\ 

in the sense that 

/ adPS irj = / (lZa)(b',b)dps irj . 

JS~D JBxB\A 

(2) On the quotient S*Xr = r\SD = Y\PSU (1, 1), we define the automorphic Patterson- 
Sullivan distributions as follows: the distribution PSf G V(S~K r ) by 

/ adPSl = / (xa)dPS irj = / TZ(xa)(b' ,b)dps irj , 

Jsx. r JSU JBxB\A 

where x is a smooth fundamental domain cutoff. To emphasize the cutoff we some- 
times write PSf (a) = ^ SIi (xo)dPSi rj . 

The following proposition is obvious from the definition, but important: 

Proposition 3.4. PSi rj is a geodesic flow invariant and Y -invariant distribution on ST) = 
D x B; and PSf r . is geodesic flow invariant on SXp. 

The geodesic flow invariance of PS irj on ST) is trivial; for PSf r . , it is also easy and results 
from the following principle: 

Lemma 3.5. Let T G T>'(T\ST)) be a Y -invariant distribution. Let a be a Y -invariant 
smooth function on ST). Then, for any a±, a<i G V(S~D) such that X^ 7 gr a i(l-(. z i &)) = a ( z i 
(i = 1,2) we have 

/ a\dT = I a2dT 

Proof. Let x be a function on Cq°(D x B) such that ^ 7 er x("f-( z i &)) = 1 ( m general, we 
choose x to be independent of b). For any such x we have 

J su aidT = f SD {j2yerX(l(z,b))} ai(z,b)T(dz,db) 

= Jsd ^ 7 er X( z , b)ai(-f(z, b))T{dz, db) 
= f S T>x(z,b)a(z,b)T(dz,db). 
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□ 

If we look at the expression 
(3.13) PS irj (a) = j Ib-b'l-^n^T^dtyT^idb'), 

and apply Otal's theorem saying that T ir . = F' ir . for some Holder function F irj , we easily 
derive: 



For any a G C°°(ST)) we have 

d 2 



\PS irj (a)\ < WF^Wi^. \\g^\b -6r Wiri ^(«)IU-(Bx B \A) 

Of course, the right side may be infinite. 

For future reference, we state a sufficient condition to obtain a non-trivial estimate: 

Proposition 3.6. Assume that \b - b']- 1 ' 2 ^ K(a) G C 2 (B x B). Then 



r\2 

\PS irj (a)\ < WF^Wl^. \\g^\b - 6'|- 1 - 2ir ^(a)|U» (BxB \ A) < oo. 

A simple example where the condition holds is where a G C 2 (ST)). In that case, there 
exist C > and K > such that: 

(3.14) \PS irj (a)\<C(l + \r,\) K \\a\\c^ 

for all j. If a G C 2 (SX r ), |PS£.(a)| < (1 + y)^ \a\ \ C 2 for all j. 

3.2. Microlocal lift and Wigner distributions. We now give a precise definition of the 
matrix elements (Op(a)<f)i rj , 0j r .) and hence of the Wigner distributions. When a is a T- 
invariant function on ST3, then in the non-Euclidean calculus ^2.21 we have 



(3.15) Op(a)cf> iri := / a(*, %^+^<*' & >T ir ,(eZ&) 



Definition 3.4. The Wigner measure of (pi r . is the distribution W$ r . on S*Xr = r\iSD 
defined by: 



a(g)W irj (dg) := (Op(a)(j) irj ,(j) ir .} L 2 {Xr) , 

sx r 

where Op(a) is given by (|3.15|) . 

To see that Wi Tj is a distribution of finite order, we note that (Op(a)(f>i rj , <fiirj) L 2 (x r ) is 
bounded by the operator norm of ||Op(a)|| and hence by a C k norm of a. In fact, Otal's 
regularity theorem shows that it is of order 2 at most. 

Wigner distributions are fundamental in the theory of quantum ergodicity. Let us recall 
the basic result: 

Theorem 3.7. |Sht IZ] Let dfi denote Haar measure on S*Xr. Then 

T77TT Y] I / adW ir — — / ddfl\ 2 -> 0, 

where AT (A) i/ie normalization factor : < A}. 
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It follows that a subsequence (Wj h ) of density one of the Wigner distributions tends to 
Liouville measure (which equals normalized Haar measure in this case). The "quantum 
unique ergodicity" problem is to know whether there exist exceptional subsequences with 
other limits. Recently, E. Lindenstrauss proved that no such exceptional sequences exist in 
the case of Hecke eigenfunctions on arithmetic surfaces jEj. 

4. Proof of Theorem 11.11 
We begin the proof witha lemma giving the explicit expression of Wi rj : 
Lemma 4.1. We have 
(4.1) (Op(a)(f) irj ,(f) irj ) L 2 {Xr) 

= 2< 1+ *"> / BxB (/ D Xa(^&)[cosh S6 ,, 6 (,)]-^)d^(.)) Tl f^^\ 

where coshsj, lj 6 2 (^) is given by The right hand side is independent of the choice of x- 

Proof. By the generalized Poisson formula and the definition of Op(a), 

(Op(a)<fe,,0 ir? .) = / (I xa(^,6)e (| " HrjKl! ' 6> e ( ^ )< *' 6 ' > ^(«) > ) T ir Ub)T ir .{db'). 



JBxB \JL> / 

We then use the following identity 

Lemma 4.2. [N] Let z G D, let 61,62 £ B and let s blyb2 (z) denote the hyperbolic distance 
from z to the geodesic 7& 1) f> 2 defined by (61,62)- Then 

, / \ 2 \ z - - 6 2 | 

COSh S blM (z) = — r-fir. 

1 61 - 6 2 |(1 - \z\ 2 ) 

Combined with ()3.10|) and ()3.8|) . we get 

e (z,b) e (z,b') = 4 [ coshsb , ]b ( z )]-2| 6 _ 6 /|-2 

Raising both sides to the power | + itj completes the proof. □ 
The next step is to analyze the integral operator 



D 



(4.2) / xa(^,6)e ( 5+-)( 2 . fe ) e (5+-)<^''')^ /( 2 ) 

= 2- {1+2ir ^ [ xa(z,b)[coshs b ^ b (z)]^ 1+2i ^\b - b'\-^ 1+2i ^dVol(z). 



D 



In this paragraph - and probably also in the rest of the paper - we sometimes drop the 
j-indices of rj, indexing the eigenfunctions by r instead. 

If we drop the factor |6 — b'\~( 1+2vr \ the right side of 14. 21 defines the operator C r : C C (D) — > 
C(B x B) by 

(4.3) C r {xa){b\b) = [ xa(z,b)[coshs b/jb (z)}- (1+2ir) dVol(z). 



D 



We now rewrite the integral in terms of coordinates z = (t, u) based on the geodesic j b ',bi 
after which we can relate C r with the operator in (jl.5|) . 
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Given a geodesic 7?/^, we work with special coordinates on D or H, adapted to 7^ as 
follows. We write z = (t,u) where t measures arclength on 7^ and u measures arclength 
on horocycles centered at b. More precisely, we denote by g(b',b) the vector on jb',b which 
is closest to the origin, and the coordinates (t, u) parametrizing z are defined by (z, b) = 
g(b', b)a t n u . For any given (&', b), the volume element of z is dVol = \dtdu\. The computation 
is most easily checked in the upper half plane, with b = 00, b' = i and g(b', b) = e = (i, 00). 
Then atn u i = e t (i + u). The area form is Setting y = e 4 , x = we* shows that the area 

form equals \dtdu\. 

We obtain 

(4.4) c r{M - /^^..--W.^)!**!. 

We further simplify as follows: 
Lemma 4.3. We have 

C r ( X a)(b,b') = [ (l + u 2 )- C * +ir) Xa{g(b,b')a t n u )\dudt\. 

JRxR 

Proof. We recall that sv,b(t,u) is the distance from the basepoint of ga t n u to the geodesic 
generated by g in the hyperbolic plane H = G/K. That distance depends only on u and has 
the value cosh sv : b(t, u) = y/1 + u 2 . 

□ 

Next, we further rewrite the operator C r in terms of the operator L r in ljl.5j) : 
Lemma 4.4. We have: 

(Oj9(a)0 ir ,0 jr ) = 2-( 1+2 ^j G L r ( X a)(g)PS ir (dg)- 
Proof. Lemma f4.3l states that 

£ r { X a){b,b') = f n L r { X a){g{b,b')a t )\dt\ 

= K(L r ( X a))(b,b'). 

Integrating against dpsi r and using the formula in Definition completes the proof. 

□ 

The next step is to apply the stationary phase method to L r ( X a). The stationary phase 
set of (|4.4|) is the geodesic jb',b from b' to b or equivalently it is the set u = in the integral 
defining L r ( X a). Stationary phase method gives the asymptotic expansion 

(4.5) L r ( X a)(g) = (-Ainr)- 1 / 2 ( £ r~ n L 2n { X a){g)) 

n>0 

where L 2n is a differential operator of order In on ST): L is the identity, the other L 2n are 
differential operators in the stable direction, that is, in the direction n u generated by the 
vector field X + . 
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If we now integrate 14.51 with respect to PSi r , and compare with Lemma l4.4[ we get an 
asymptotic expansion, 

(4.6) (Op(a)0 ir ,0 ir ) = 2-( 1 + 2 -)(-4mr)- 1 / 2 (Vr-" / L 2n { X a)(g)PS ir {dg)) 

n>o Jsii 

Because the distribution on the left-hand side, e^ +ir) ^ b) e^ +ir){z ' b,) dVol(z)T ir (db)T ir (db'), 
is T-invariant (as a distribution in the triple (b, b', z)), each of the distributions obtained in 
the expansion, 



/- / L 2n (f)(g)PS ir (dg), 
Jsr> 

is T-invariant. In application of the principle IH.5| the functional 

fl H L 2n (xa)(g)PSir(dg) 
Jg 

defines a distribution on T\G, and the definition does not depend on the choice of \. The 
first term (n = 0) is precisely the Patterson- Sullivan distribution PSi r as defined in the 
quotient 5X r . 

4.1. Completion of Proof of Theorem 11.11 We now turn to the relation between Wj r 
and PSi r . It follows from the stationary phase asymptotics above, (|4.5|) . that 

/ a( 9 )^ J (d S ) = 2-('*>fr 7 »/ WM ,^t tO(r;''- tlf ) 

where K was defined in KS. 141 If we choose N > K then the remainder term goes to zero. 
Since L = Id, the operator L r N ^ = ^2n=o r ~ n ^ 2n can ^ e inverted up to 0(r~ N ~ 1 ), that is, 
one can find differential operators Mr N ^ = Y2n=o r ~ n M 2n (with Mo = Id) and Bi N ^ such that 

L r N )M^ = Id + r- N - l R r N \ 

We thus get 



sx 



MiPa( 9 )W iri (d9) = / sd + 0(rj 



JV-1+A"\ 



The second line is a consequence of 13.51 Since we know, from standard estimates on pseudo- 
differential operators, that the Wigner measures are uniformly bounded in (C k )* for some 
k, we have 

M^a(g)W irj (dg) = j a(g)W irj (dg) + 0{rj l ). 
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This shows that 



r l+2ir) I a(g) PS : r + d9 l = j a{g)W ir Ug) + 0(rf) 



5x r {-Amrj) 1 / 2 Jsx 



and finishes the proof of Theorem 11.11 



5. Integral operators and eigenfunctions 

In this section, we give further results on the operators L r (11.5)1 and Ips ir (jl-7j) which 
will be needed in the proof of Theorem 11.31 With no extra work, we treat general integral 
operators of the form 



(5.1) !M(u) := / cr(gn u )fi(dg), 

Jr\G 

where a G C°°(T\G) is an automorphic form and where d\x is an invariant distribution for 
the geodesic flow on T\G. In addition to \i = PS ir the case where /i is a periodic orbit 
measure is also important in this article. In this case, we write 1^ = J 7 with/ 7 (<r)(n u ) = 
J i j a(a~ l an u )da. Here, a 7 G G is an element conjugating 7 G T to an element of A. This 
expression arose in the trace formulae of [Z] and will arise in ^1 The similarity of these two 
kinds of integral operators may be seen as one of the deus ex machina behind Theorem 11.31 

5.1. The integral operator J^. We can view 1^ as an integral operator from C°°(T\G) — > 
C°°(N) ~ C°°(R) The following lemma shows that when a is a joint eigenfunction of the 
Casimir operator and of W, then 7 M (cr) solves an ordinary differential equation in u. When 
a is a (r, m)-eigenfunction in the complementary or principal unitary series, the equation is 

(5.2) (u 2 + 1)±L + (2u -»)| + (- + r 2 )/ = 

We denote by F T>m (^=|) the even solution of (|5.2|) which equals 1 at u — 0, and by G Ti7n {^i) 
the odd solution whose derivative equals 1 at u = 0. In the holomorphic discrete series, and 
when a is the lowest weight vector, the analogous equation is the first order equation 

A basis for its solutions is given by f(u) = (—i)~ m l 2 {u + i)~ m l 2 . There are similar equations 
for higher weights and for the anti-holomorphic discrete series, but for simplicity we only 
discuss the lowest weight case. 

Proposition 5.1. Let fi be a geodesic flow invariant distribution on T\G. 

• Let a be a (r, m)- eigenfunction in the principal or complementary series. Then 
I^(o~)(u) is a solution of 115.0)) . Hence, 



F ( u — % \ / / / it — z 

I^(cr)(u) = ( J adjj)F T)m \^Z%J + ( J X + adji)G T ^ m \^ 



-2i 



where F,G are the fundamental solutions of defined in [2] (2.3) (see \5. 7| ) for 
formulae in terms of hypergeometric functions)). 
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• Let a be a (T,m)-eigenfunction in the discrete holomorphic or anti-holomorphic se- 
ries. For simplicity, assume a = ip m (the lowest weight vector in T>^). Then: 

= ( / ^)(-o~ m/ >+*r m/2 - 

Jr\G 

Proof. In the case of 7 M = J 7 , the proof is given in [Z], Proposition 2.3 and Corollary 2.4. 
We briefly verify that the same proof works for any invariant distribution. 

First assume cr is a (r, m)-eigenfunction in the continuous series. Since 4f2 = H 2 + 4X+ — 
2H - AX + W we find that 

(5.4) (a^- 2 - Aim^- + 4(1 + r 2 )) I,{<j){u) = - J {(H 2 - 2H)a) (gn u )^dg). 

We write Ho~(g) as 2-^ t _ a(gat)- Using that n u at = atn ue -2t and that \i is an A- invariant 
distribution, we find that f Ha(gat)^(dg) = —2u-^I fl (a)(u). A similar calculation replaces 
H 2 by the square of this operator. The final equation is as stated above. We then evaluate 
J M (cr) and its first derivative at u = to obtain the expression in terms of F, G. 

In the discrete holomorphic series, we use that E_a = to get 2iX + a = (H — m)o~. This 
leads to equation ()5.3|) and to the solution given above. 

□ 

5.2. The integral J R (1 + u 2 )~ s I il (o~)(u)du. In Theorems 11.21 fOl and elsewhere, we will 
need explicit formulae for the integrals 

(5.5) / {l+u 2 )- s I^{a){u)du 

Jr 

We assemble the results here for future reference. 

In view of Proposition 15 .11 we need explicit formulae for the integral of (1 + u 2 )~ s against 
the functions F r , m (^=f) , G T ,m (z§{) , and (-i)- m ' 2 {u + i)" m / 2 . In fact, by Proposition O 
and 12.51 it will suffice for Theorems 11.21 and 11.31 to have explicit formulae just for F t0 and 

(-i)- m ' 2 (u + i)- m i 2 . 

We introduce the following notation: 





:= 2 f™(u 2 + l)- s du, 


/4 fc ( s ) 


:= J R (u 2 + l)- s F(\ + 2 -^,\- 


< /4,m( s ) 


: = r oo ( M 2 + l)-F r , m (^|)^, 




:= J R (u 2 + l)- s G T , m (^f) du, 




:= J R (u + i)~ m/2 (u 2 + l) s du, 



It is clear that the integrals defining /x (s) and /i„(s) converge absolutely for s > \ 
and 9ft(2s — y) < —1, respectively. We now show: 
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Proposition 5.2. The integral defining [i c irk {s) converges absolutely for — 23?s — l/2 + 9farfc < 
— 1, and in this region we have: 



/oo 
(M + 1) -2* S -l/2 + ^ 
-oo 



for some constant C (independent of s,rk)- 

Proof. Indeed, as in [Zj (Proposition 2.7), the differential equation ()5.2|) is equivalent, by a 
change of variables, to a hypergeometric equation, and a short calculation shows that 



(5.7) 



rp ( u-i \ _ rp(l I 2ir k 1 _ 2ir k 1 
r ir k ,0\_2i) r \Z> 4 '4 4 '2 



,-« 2 ), 



r< / ti— i \ _ 17/3 i 2ir k 3 2ir k 3 ,,2 



' 2' 



Classical estimates on hypergeometric functions (see also Z , p. 50) show that there exists 
C > (independent of r k ) such that 



p/l I 2irfc 1 _ 2ir k 1 
^ V4 ~r 4 ' 4 4 ' 2> 



-It 



(5.8) ^ <C(1 + |«|) 



-l/2+5Rir fe 



p/3 , 2irfc 3 _ 2irfc 3 _ 2 
1 \4 T 414 4 ' 2 ' 



These estimates follow immediately from the connection formulae for hypergeometric func- 
tions: 

F(a,b,c,z) = r§$E$(-z)- a F(a,l + a- c, 1 + a - b; z' 1 ) 

+ W0)^ bF ^ l + b-c,l + b-a; z" 1 ). 
Since F(0) = 1, we obtain that (as \u\ — > oo) 

p/l ,2^ 1 _ 2irfc 1 2\ ^ r(|)r(-ir fc ) , , ( l + ^ rfe ) 

1 U T 4 M 4 ' 2' 14 / r(i-^) 2 11 

(5.9) 

The asymptotics ()2.1())1 of the T function on vertical lines shows that the ratios of T func- 
tions are uniformly bounded in rv The decay rate |u| _ ^2~ M ' fe ) is sufficient for the absolute 
convergence of the integral in (|5.6|) as long as 9ft(| — ir&) > 0, i.e. if ir^ is not the parameter 
of the trivial representation. 

□ 



Although we will not need them, we note that the estimates for G are similar. Each of 
the above functions admits a meromorphic continuation to C. Since we will not need the 
results for general /i£ m (s), ^ c °^{s) we omit them in the following. 

Proposition 5.3. We have: 
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AioOO = r % ( ^ } (Ke S >|) 

, /C r _x r(|)r(s-i+^)r(s-i-^) 

^ir k \ b ) — r(s) 2 

d , x (-i) m/2 7r2 2s+2 - m / 2 r(-2s+f ) 



(3fte s > 0) 



* m ^y _( 2s+ i-Ht)r(-s)r(-s+f ) (^(2s 2 )< 1). 

The proof is given in [Zj (see pages 50-52). 

6. Proof of Theorem 11.21 
The key objects in the proof of Theorem 11.21 are the closely related integrals 
f I r (a) = J R (l + n 2 )-(l+-)p^ r ((a X )")^, 

(6.1) 

where PS iTj G V'(B),PSf r . G V{T\G) are defined in Definition IT31 and where f u (g) = 
f(gn u ). Note that PSf r (o~ u ) = I PSir (cr)(u) in the notation of §5] It takes some work to 
prove that each integral is well-defined. In Lemma lb. II it is proved that the two integrals are 
well-defined and equal for a G C°°(T\G). 

Theorem 11.11 equates the Wigner distribution with the distribution a — > j su {L T:j xo~)dPSj 
for a G C°°(r\G). In Proposition lb .41 we show that this functional also equals I r (cr) = 
The explicit formulae for the Wigner distributions in terms of the Patterson-Sullivan 
distributions follow from the identification with i^(cr), which can be explicitly evaluated 
using the results of £0 

6.1. Convergence and equality of the integrals. In the following, we recall that ^tir^ = 
in the unitary principal series but is positive in the complementary series. 

Proposition 6.1. We have: 

(1) Integral converges for all o G C°°(T\G) ifVi r is in the unitary principal series. 

(2) I r {a) = lf{cr) for all a G C°°(T\G). 

6.1.1. Proof of (1). We give a representation theoretic proof that 
(6.2) f (l + u 2 r {1 > +ir) I T PsM(u)du 

converges absolutely. We make no attempt at a sharp estimate but only one sufficient for 
the purposes of this paper. 

Lemma 6.2. Let PSi r be the Patterson- Sullivan distribution corresponding to <f> r . Then: 
(i) I PSir .(4> irk )(u) < C (1 + |r fc |) 4 (l + \ rj \) K (l + luD-i/a+Kto). 

(it) I PSiT .&m) <C(1 + |m|) 4 (l + \ rj \) K (l + \u\y m , 
where K is the same as in \3.1$ . 
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Proof, (i) By Propositions 15. II and 12. 5^ 

Ips ir . OirJ (u) = J nG irfc (gn u )PS irj (dg) 

(6-3) 

= {Ir\G^ k PSi rj (dg))F irk>0 [^) 

By (jH.14j) . there exists .ff so that 



m rk ,PS trj )\<C(l + \r,\) K (l + \r k \) 



Here, we used a crude estimate ||0irj|c 2 < C(l + l r fc|) 4 (in fact, rl/logr^ is true, but it is 
not necessary for our argument). We combine with the estimates in Proposition 15.21 fcf. 15. 9\i 
on the hypergeometric factor to obtain the estimate stated in (i). 

(ii) We now have 

Ips iTj (tpm)(u) = f r \ G ^ m (gn u )PS irj (dg) 

(6.4) 

= (J^ G ^ m PS irj (dg))(u + i)- m / 2 . 
To complete the proof we note that \(u + i)~ m / 2 \ < C(l + M)~ 1/2 and that (by (ETT3II ). 



\(^ m ,PS irj )\ <C(l + |r,-|) K (l + |m|) 4 



□ 



Given a co-compact discrete group T C SX(2, R) we denote by r = 3(iro) the Casimir 
parameter corresponding to the lowest non-zero eigenvalue of A, i.e. the complementary 
series representation closest to the trivial representation. 

Lemma 6.3. If c E C°°(r\G) has no component in the trivial representation, we have: 



I PSir (a)(u)<C(l + \r\) K (l + \u\) 



-l/2+ro 



Proof. Since PSi T is geodesic flow and time reversal invariant, we may write by Proposition 



(6.5) W*)(«) = X) 
It follows by Lemma (6.21 that 



IW*)(«)| < C(l + |r|)^x{V (1 + M) 



(6.6) 



(1+ 



EJi + H) 4 l^*yl (i + H)- m ) 



By Proposition 12.41 



, — ir ,■ ) 



< ||X_7> ir »lb, ly^^r 

X^Wj — m/ 



It follows that for any M > there exists a constant Cm so that 



(6.7) 



1 ) 



<C M (l + |r,|)- M , 



<C M (l + |m|) 



-Af 
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Indeed, 

V ir {cr) = Vir,m<f>ir,m, with |er inm | < C M (1 + \rj\ + |m|)~ M , 

mGZ 

hence 

\\X_V irj (a)\\ L 2 < C M ^2(l + \rj\ + \m\)- M (l + \ rj \ + |m|), 

m 

where we bound | |X_0 ir . jm | | L 2 < C(l + |rj| + |m|). 
Similarly, 

oo 

(P^)(a) = CT m , m +2n^m,m+2n, With \a m ^ m+2 n\ < C M (l + \m\ + \n\)~ M , 
n=0 

hence 

\\X_Vi(a)\\ L i < C M ^(1 + |m| + |n|)- M (l + |m| + |n|). 

m 

By f|6.7|) and Lemma 1)6.20 . the sum ()6.5|) converges absolutely and the decay estimates in 
u sum up to the stated rate. □ 

Completion of proof of Proposition 16.11 (1): It follows from Lemma f6. 31 that 

|/ r r (<r)| < J R {l + u 2 )-^ +m ^\\PSi{o- u )\du 

< C(l + |r|) x / R |(l + u 2 )-^+^)|(i + \ u \)-y^du. 



(6.8) 



Since Vi T is in the unitary principal series, Jfor = and so |(1 + M 2 )~^+ ,r )| = (1 +-u 2 )~2 and 
since — 1/2 + r > it follows that the last integral in (|6.8|) converges absolutely. 
We now move on to the assertion (2) of Proposition 16.11 

6.1.2. Proof of (2). By Proposition 13.51 we have 

(6.9) / a(gn u )x(gn u )dPS ir (g) = / cr(gn u )x(g)dPS ir (g). 

Jg Jg 

Indeed, x(sO an d o\ u {g) '■= x(.9 n u) are both smooth fundamental cutoffs, so both sides equal 
f r ^ G adPSf r . Integrating against J R (1 + u 2 )~^2 +ir ) completes the proof. 

6.2. Continuity of PSi r . As mentioned above, the Wigner distribution equals the func- 
tional a — > Jg^Lr-x^dPSir-.. To prove that this also equals I r (a) = l£ (a) we need the 
following continuity result for the functional PSi r . 

Lemma 6.4. PSi r G T>'(ST)) has the following continuity property, 

PS ir (L r (xcr)) = [ {l+u 2 )-^PS ir {ax) u )du, 
Jr. 

where f u (g) = f(gn u ). 
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Proof. By Definition E31 

PS t . r ((ax) u ) = ps ir (K(a X ) u )) 

= faxvUru ^ ^ b')a t n u )dt} d ^00& 
We first note that for all u, TZ(xa) u E C™(B x B' \ A) since (x^T e C C °°(5D). It follows 
that pSi r (7?.((Tx)") is well-defined and smooth in u. 
The continuity statement is equivalent to 

(6.10) ps ir {K{L r {xo)) = J R (l+u 2 )~^ps ir (n(ax) u )(b,b'))du, 
or equivalent ly 

(6.11) ps ir ( [ (1 + u 2 )-^ +ir) K{a X T{b, b'))du = [ (1 + « 2 )" ( ^ <r W(^W)(^ b')du. 

Jn Jn 

We must again check that both sides of (|6.1U|) are well-defined. Clearly, 1Z(L r (x°~)) is 
well-defined because x nas compact support. The problem is to prove that the left-hand 
side is well-defined, since that psi r is only known to be a bounded linear functional on 
| b — b'\ 1+2ir C 2 (B x B' \ A) (cf. Proposition 13. 6|) . We therefore have to verify that 

/ (1 + u 2 )-^ +ir) K(a X ) u (b, b'))du 6 \b - b'\ 1+2ir C 2 (B x B' \ A). 
Jn 

By Lemma [4.31 and ()4.3j) . we have 
(6.12) 

(l + « 2 )-d+ <r )^( ( Txr(6,6 , )) d «= |&-&'| (1+2ir) / (x^)(^&)e ( ^ +ir)<z ' f>> e { 5 +ir )^'>d\/ /(z), 



and therefore the condition to be satisfied is that 

(6.13) [ (x<?)(z,b)e {1 i +ir){z > b) e {1 i +ir){z > b,) dVol(z) G C°°(5 x £). 

Jd 

This is clear due to the compact support of x° lYy z i which is independent of (b, b'). 
We may then rewrite ()6.11|) as: 

((/ D (Xff)(^,6)e(5+^)^)e(H^)(^')dKo/(z)),dT, r ® rfT ir ) 

(6.14) 

= f D X ({o-(z, b)e^ +ir ^e^ +ir ^ z ^), dT ir ® dT ir )dVol(z), 

where (, ) denotes the pairing of test functions and distributions on B x B. I.e. we need to 
check that we can pass dTi r ® dTi r under the dVol(z) integral sign. 

By Otal's regularity theorem (see Theorem 13. 2|) . dTi r = F^(b)db where Fi T is a continuous 
27r periodic function in the sense that Fi r {Q + 27r) — Fi r {9) = C r . Integration by parts then 
gives 

(9,T ir )= [ g(b)T ir (db) = - [ g'(b)F tr (b)db + g(0)(F ir (2n) - F ir (0)). 

J B JB 

Applying this in each of the (b, b') variables to the pairing on B x B in (|6.14j) produces four 
terms of which three involve the boundary term (Fj r (27r) — F ir (0)) and the fourth is 

{ / d 2 (( X a)(z,6) e ^ + ^^e(l +ir )^ 6 '>dF O K^))}F ir (6)F ir (60^®^'- 
BxB Jd Co x oo \ J 
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By applying Fubini's theorem to the fourth term, we obtain: 

IbxbUd a&W b)e^^ b h^^dVol(z^ }F ir (b)F ir (b')db ® db' 

(6.15) 

= Id xUb.b a&F b)eG++M>e$+*M) T ir {b)T ir {b')dbdb'}dVol{z), . 



Fubini's theorem applies in a similar way to the other terms. We then transfer the b deriva- 
tives back to dT ir and obtain (jfi.!4j) . 

□ 

As a corollary of Proposition 16.11 we obtain the following explicit formula: 

Corollary 6.5. We have: 
(6.16) 

+ E^^jl/^W^Wg + ir). 
All integrals and series converge absolutely. 

Proof. In fact, by Lemma 16. HI we may interchange the order of summation in (j6.5j) and 
integration in (J6.2j) . Using (J6.3j) and ()5.6j) . we have 

/ (1 + u 2 )-^I T PS%r {<p ir] ){u)du = ( / ^PSUdg))^ + ir), 
Jn Jr\G 1 

and thus obtain the first series of the stated formula. Using ()6.4)1 and 1)5.6)1 . we have 
(1 + u 2 )-^I r PS J^ m )(u)du = ( / ^ m PS ir {dg))it{\ + ir), 



r Jt\g 2 

and obtain the second series. □ 

6.3. Completion of Proof of Theorem ll.3L To complete the proof it suffices to explicitly 
evaluate I^{o~) for the generating automorphic forms. 

Lemma 6.6. In the special cases when a = <pir k , X + <ft irk or ip rn , we have the explicit formulae: 

(1) In the case a = (j) iT . k , (Op(0 ir j0 irj , (f) irj ) = A4 fe (§ + ir j)(f 4>ir h ^P 'Sf r .(g)) . 

(2) For a = X + <p iT . k , (Op(X + (fr irk )(f) irj , (j) ir .) = for all j . 

(3) For a = ip m , (Op(ip m )(j) ir . , <fi irj ) = ^(\ + ir s ) ( / ip m dPSf r . (g)) . 

Here, the expressions n\ r (~ + irj), + i r j) are defined in \5.b}) . 

Proof. All statements follow from the combination of Proposition 15.11 and Proposition 15.31 
The case a = X + (p irk follows from Proposition 12.51 and the fact that the Patterson-Sullivan 
distributions are invariant under time- reversal (cf. Proposition 13.3)) . 

□ 

We note that these explicit formula give a new proof of Theorem 11.11 

Corollary 6.7. When a is a joint (Q,W)-eigenj 'unction, we find again that {Op(a)(pi r ,<f)i r ) 
is asymptotically the same as r -1 / 2 J adPS ir . 
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Proof. By definition, F Tjm (l/2) = 1 whereas G T>m (l/2) = 0,G"(l/2) = —2i. The stationary 
phase method then shows that 



whereas 



Here, we use the estimates in ()6.2|) . which are valid in all weights. 

7. Dynamical zeta-functions: Thermodynamic formalism 



□ 



In this part, we relate the Patterson- Sullivan distributions to the Ruelle resonances intro- 
duced in |R"u87j. 

Ruelle's aim was to study the Fourier transform of the correlation function, 

p FG (t) = J F(x)G{g t x)duo{x) - J Fduo J Gdu, 

(t > 0), in the very general context of an Axiom A flow (g l ) (for instance, when u is the 
measure of maximal entropy). He showed, for smooth enough functions F and G, that the 
Fourier transform pp,G has a meromorphic extension to a half-plane of the form {3?es > h— e}, 
strictly beyond its half-plane of absolute convergence {3?e s > h} (where h represents, in a 
general context, the topological entropy). He used the so-called "thermodynamic formalism" 
and showed that the poles of p(s) occured precisely for certain values s, linked with the 
existence of certain distributions obeying specific transformation rules. 

In the case of the geodesic flow on a compact surface of constant curvature —1, and for 
C l functions F, G on T\G, the Fourier transform p is an analytic function in the half-plane 
{s, Kes > 1} and has a meromorphic extension to {9ftes > 0} with poles at values of s = h+ir 
for which there exists a distribution Ci r on S*Xr satisfying: 

• g l .e ir = e- ( ^ +ir)t e ir 

• 6i r is invariant under the stable horocycle flow. 

In the case of constant negative curvature, it follows that e, r is given by: 

(7.1) e ir (F) = [ F(z,b)e {1 z +ir){z ' b) T ir (db)dVol(z) = [ Op{F)<t> ir (z)dVol{z) 



where T ir is the boundary values of an eigenfunction <p ir of A of eigenvalue \+r 2 (see equation 
(13. 5|) . and [Z2] ) Hence the poles of p, i.e. the Ruelle resonances, occur at s n = 1/2 + ir n . If 
the eigenvalue is simple, the residue of p a ^ at s n is given, up to multiplicative normalization, 
by 

R rn (F,G) = ( [ F(z,b)e il/2+irn){z ' b) T irn (db)dVol(z)) ( [ G o L (z, b)e {1/2+irn){z > h) T irn (db)dVol(z 



Op(F)(/> irn dVol(z) / Op(G o L)<t> irn dVol{z) ) = e irn (F)e irn (G o t) 
x r / \</x r 
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where i denotes time reversal. To see this, we observe that the residue R rn (F, G) is bilinear 
in F and G, and its definition implies that it satisfies the identities 

R rn (F o g\ G) = R rn (F, G o g~*) = e -<Wr»)« Rr n (F, G), 

and 

R rn (Foh u + ,G) = R rn (F, G o K a _) = R Tn (F, G), 

where h+ denotes the stable horocyclic flow and h u _ the unstable one. It follows that it must 
equal ei rn (F)ei rn (G o i) is the eigenvalue 1/4 + r 2 is simple. If the eigenvalue is not simple, 
the expression becomes more complicated, as one has to form a linear combination of the 
functional associated to the various eigenfunctions. 

We now prove prove Theorem 11.31 for Z^- In this section, the Patterson-Sullivan distribu- 
tions are normalized so that J UPSi Tn = 1. 

Remark 7.1. The residue is obviously not the same as in Ruelle's paper: in our case the 
residue must define a geodesic flow invariant functional, whereas in Ruelle's paper it has to 
be invariant under the horocycle flow. 

7.1. Markovian coding of the boundary. The proof given here relates the function Zi 
to the determinant of certain operators, called transfer operators. To define them, we need 
to recall from jSej the construction of Markov sections, using the Bowen-Series coding of the 
action of T on the boundary B. Series used this construction to study Poincare series. We 
apply it to the somewhat different objects Z 2 . For this application, we need some further 
discussion of Markov coding which we could not find in the literature. 

If we want to study the action of V on the boundary, and the existence of conformally 
invariant distributions - by this, we mean the property 13 .41 - it is enough to consider a set of 
generators of V. In fact, it is even enough to work with a single, well chosen transformation 
of the boundary: roughly speaking, this transformation F^ will be defined by cutting the 
boundary B into a finite number of closed intervals Jj, and will act on each Jj by the action 
of one of the chosen generators of V. 

We will require the map F^ : J = UJj — > J to have the following properties: 

(i) F^ is analytic, expanding (or at least, some power of F^ r ' is expanding). 

(ii) The Js form a Markov partition for F^ r \ This means that F^ r ' sends the boundary 
of J to itself. 

(iii) The natural map J = U Jj — > B gives a bijection between periodic points of F^ and 
points at infinity of closed geodesies, except for the closed geodesies ending on the boundary 
of an interval Jj, that have exactly two preimages. If F^ r ' x = x, and 7 is the closed geodesic 
corresponding to x, then \(F^ )'x\ = e Ll . 

We recall briefly the construction of F^ proposed by Series jHij, when T is cocompact: she 
chooses a symmetric generating set for T, called T . She then defines a notion of "admissible 
representation" of an element g G T as a word g = g\...g n with gj e r , such that 

- an admissible word is a shortest representation of g in the alphabet r . 

- every g G To has a unique admissible representation. 

Without going into details, admissible words are shortest word-representations; and be- 
sides, whenever there is a choice of several such representations, one chooses the one that 
"turns the furthest possible to the right" in the Cayley graph of T with respect to To (seen 
subset of the hyperbolic plane). 
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Let us denote Trp the set of finite admissible words; the notation is borrowed from [Hi] 
but we are adding an r to specify that we are choosing representations that turn the most 
possible right in the Cayley graph - the same convention Series used in her paper. Replacing 
"right" by "left" one would obtain another notion of admissible words, and we denote Ey 

the set of left-admissible words. Note that E^P = E^ . Now define E^, the set of infinite 
right-admissible words, as 

E« = {{ 9j ) e if .//,,..//,., e Z ( ;\Vj,k> 0}. 

Series shows in jSi] that the map 

(7.2) Sj. r) — > H, g l ...g n h^ gi ...g n .O 

can be extended to a continuous map j( r > : E^ — > B. She denotes I^(gj) the set of points 
in B that have a representation by a sequence in E^ starting with the generator gj. The 
boundary B is thus divided into a finite number of closed intervals, with disjoint interiors. 
One can define a map F( r ' that acts on E' r ' by deleting the first symbol and shifting the 
sequence to the left; seen as a map on B, it acts as g~ l on each interval I^(gj). Actually, the 
map is defined on 1^ := Ul^(gj); when working on B one should always remember 
that its definition is ambivalent on boundary points. The partition B = Ul^(gj) is not 
exactly a Markov partition for the action of F^ r ': there is no reason that boundary points 
should be sent to boundary points. But, by construction, the images of these boundary 
points under iteration of F^) form a finite set. Cutting the intervals I^ r \gj) at these points, 
one can refine the partition B = Ul^(gj) into a new finite partition B = UJj that is now 
Markov. This way we obtain a tranformation F^ satisfying all the conditions (i), (ii), (iii). 
An element of B may be coded either by a word in £W, as we have already seen, or by an 
element of the subshift of finite type 

A (r) = {(i k ) k > ,F^(J ik )nint(J ik+1 ) for all k > 0}. 
Both codings are bijective except on a countable set (in fact the coding map is at most 2 to 
1). " 

To make the link with the geodesic flow, we now extend the expanding transformation 

F (r) 

to an invertible transformation F of a subset of B x B; in terms of symbolic dynamics, 
we want to work with two-sided subshifts. We consider 

E« = {(&■) G \f.g,...g,., G E^ , Vj, k > 0} 

and 

£(0 = {(&) e if •//, /,...//, e 4 r) Vj,* < °>- 

Formally, elements of E(j) are inverses of elements of E^). By the same consideration as 
before, we have a coding map j"' : E^ — > B or equivalently jm : E/n — > B; we denote 
I(i){gj) C B the interval formed by points whose coding in Em ends with gj. This gives a 
partition B = \Jln\(gj) and a map Fn) on B, that corresponds to the shift to the right on 
Em. If we want we can refine the partition B = Ul^(gj) into a Markov partition B = UKj 
and code the dynamics by a one-sided subshift of finite type A^. 
Let us now introduce the two-sided subshift E, 

E = {(g^ G r%,gj...g j+k G Yiy,Vj,k}\ 
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E is naturally isomorphic to a subset of Em x Ev), and thus there is a coding map from E 
to a subset X C B x B: 

j:E — > X, 

ifao.* (r) ) = toofao).i (p V r) )) 

The shift to the left on E gives an invertible transformation F on X; note, as above, that 
F is actually well defined on a subset of Ulm(gj) x Ll/M (gr^) and is defined ambivalently 
at certain points of X. If y is in I^ r \gj) then F(x,y) = {g7 x x, gj x y), or in other words 
F(x,y) = (Gyyx, F^ r 'y), where Gyy is the inverse branch of F^ taking values in I^(gj). 

The partition of X into Ul^(gj) x U/ (r )(g fc ) is not a Markov partition for the action of F, 
but X = U(Kj x Jj) is Markov. If we want to work with a subshift of finite type, we should 
consider the set 

A = {(m k ,n k )kez/F(K mk x J n J n int(K mk+1 x J„ fe+1 ) 7^ for all fc G Z}. 

We can identify X with a subset of the unit tangent bundle T\G that forms a transversal 
section for the action of the geodesic flow: we observe that, for each (x, y) in X, the geodesic 
7 Xj3/ contains a unique vector, denoted v x>y , which is in the stable manifold of a vector based 
at 0. To recover the action of the geodesic flow on the whole tangent bundle, one needs to 
add a time parameter measuring the time it takes to go from (x, y) to F(x, y). Because j x>y 
and jF(x,y) represent the same geodesic in the quotient T\G, there exists a unique r(x, y) G R 
such that g T v X)V = vf(x, v )- More precisely, the function r is defined without any ambiguity on 

o o 

\-\Kj x Jj and can then be extended to a continuous function on UKj x Jj. By construction, 
the function r is actually locally constant on stable manifolds; i.e., it depends only on the 
variable y. It is analytic on each rectangle Kj x Jj. We see r as a return time from the 
section X to itself; note however that r may change sign: we are not exactly in the usual 
situation of a "first return time". Nevertheless, when y is a periodic point of period n of 
FW, and 7 is the corresponding closed geodesic, we have Y^k=o t(F^ y) = J 7 > 0. 
We have a surjection, almost one-to-one, from 

X T := {((x,y),t) e X xR,t e [0,t(x)}} 

to the unit tangent bundle, defined by saying that the image of ((x, y), s) is g s v x>y , the image 
of v x>y under time s of the geodesic flow. This surjection is not one-to-one on boundary 
points; by definition of r, ((x,y),r(y)) has the same image as (F(x, y), 0). On X T the 
geodesic flow corresponds the translation of the parameter t. 

7.2. Transfer operators. Let a be an analytic function on S'Xp. Let A be the real-analytic 
function on X, defined by 

r-r{z 2 ) 

A(zi,z 2 ) = / a(zi,z 2 ,t)dt 
Jo 

if (zi,z 2 ) G X C B x B. In other words, A is the Radon transform A = lZ(x a ) defined in 
EH and x is the cut-off function x{( x >y), s ) '■= 1^2)1^ (24 )l(o,r(y))(>)- If z = (zx,z 2 ) 
is a periodic point of period n for F, and if 7 is the corresponding closed geodesic, then 

71— 1 „ 

S n A(z) :=J2MF k z) = a. 

k=0 ' 
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Following |Rugh92| , we introduce a family of transfer operators L S}Z (s, z G C), acting on 
xiyji,^) which are analytic with respect to u) 2 and analytic functional (distributions) with 
respect to We refer to |Rugh92| for the precise definition of the Banach space to work 
with. The transfer operator L sz acts on x by 



/ 



(if) being an arbitrary analytic test function.) The sum runs over all k_i,ko such that the 
rectangle K k x x J ko contains a point (o^o^) with F(co>!, u^) = (zi,z 2 ). Here the notation 
g fco means the element of T such that F(uj x ,uj 2 ) = (g k ^u l , g^uj 2 ) if (co>i,co> 2 ) e iffc^ x Jfc - 

7.3. Determinants and zeta functions. Apart from the introduction of the weight A, 
our transfer operator also differs from Rugh's by the terms g^i^Y 3 ^^ 1 ^ and g^ 1 ' (u-i) 8 ' 2 . All 
his arguments can be adapted with obvious modifications to this situation. Let us introduce 
the notations T\(uS) = —loggZ'(uJi) and t 2 (oj) = loggj^c^) if uj — (lo\,lo 2 ) G X with 
^2 G Jk - 

In paragraph 4.4 of his paper, Rugh shows that L sz is a nuclear (trace class) operator. 
One can take the determinant of / — L s z : 

d(s, z) := det(I - L s>z ) = J](l - A«) mW 

where the product runs over the spectrum of L SsZ: and = m!t} z denotes the multiplicity 
of AW. The eigenvalues do not necessarily depend analytically on (s,z), as the multiplicity 
may vary; the determinant d(s, z), however, is an analytic function of (s, z), as follows from 
the next paragraph: 

For arbitrary (s , ^o), consider, for every i, a neighbourhood Vi of Xslz , such that the V^s 
are all pairwise disjoint. Let P\ z be the spectral projector on Vi for the operator L S;Z : P\ z 

depends analytically on (s,z), in a neighbourhood of (sq,zq). Call B\ z = L SjZ P$: they are 

operators of rank m^l ZQJ depending analytically on (s, z) in a neighbourhood of (s , z ). By 
definition the spectrum of B^\s,z) is contained in V{. Of course, 



det(l-B®„„) = (1-A« 



Sq,zqJ \^ ' K s ,z ) 

Thus one can write, in the neighourhood of (s , z ), 
(7.3) d{s,z) = Y[det{l-B {l) 



(i) 

! 0. 2 



s,z) 



This shows that the determinant d(s, z) is an entire function, and has zeros exactly when 
L S)Z has the eigenvalue 1. 

Rugh shows that the following remarkable identities hold: for all n, the trace of L n sz is 

zS n A(z)-sS n Tl{z) + (l-s/2)S n T 2 {z) 

(7-4) Tr(Ll z )= E - 



r ; ; \det(DF-(z) - 1) 



It follows that 
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Proposition 7.1. 

/ i e zS n A(z)-sS n r 1 (z)+(l-s/2)S n T 2 (z) 

(7.5) d(s, z) := det(I - L s , z ) = exp \- \det(DF n (z) — 1)| 

In particular, the function 

— ( 0) - - 1 S n A(z) 

d I s ' ) - Z_> ndet\DF n - 11 

has poles exactly when 1 is in the spectrum of L Sj o- 

Because periodic points of F correspond to closed geodesies, we can express 17.51 in terms 
of periodic geodesies. If F n z = z and 7 is the corresponding closed geodesies, we have 
S n A(z) = fa, S n Ti(z) = S n T 2 (z) = Lry. Thus, d(s, z) is more or less the same as 

/ -1 p(z J a— (s— 1)L 7 ) \ 

-p -e e sinh(pV2)P = n n a - ^-—^). 

\ 7 P>1 1 V 7/ n / 7 m,neN 

A "prime" following a sum or a product means we are summing over primitive closed orbits. 
Otherwise, we sum or take the product over all closed geodesies. 

The previous formula, however, is not exactly true, because certain periodic geodesies 
correspond to several different periodic orbits of F; namely, those going through the boundary 
of X (there are a finite number of them). The precise expression of d(s, z) in terms of closed 
geodesies given in |Rugh96| , or jM^: 

(7.7) d(s,z) = Y[ Y[ (i_ e *J>-(*+m+«K<-y) .p( Sj< z) 

7 m,n£N 

where the correction term is 

n c n m , n > a-e- (s+m+n)l{c)+zLa ) . 



(7.8) P(s,z) 



n c/ n m , n >o^-e- {s+m+n)l{c ' )+zL ' a y 



the products run over a finite number of periodic orbits that are counted several times in 
the Markov coding. The correction factor on the right is analytic and non-vanishing in 
{3?e s > 0}, thus the zeros of the two functions (|7.5|) and (|7.6|) are the same there. 

Remark 7.2. In the half-plane {9fte s < 0}, the correction factor P(s,z) is more difficult to 
analyze because it seems that its singularities could depend on the choice of the Markov 
section X. It was, however shown in |Rugh96| that the apparent singularities of (|7.6|) . 
appearing from the identity (|7.7|) . are removable. 

Remark 7.3. For z = (which is the case treated by Rugh) we obtain the relation 
(7.9) d(s,0)= ]JCs(s + n) .P(s,0) 

where (s is the Selberg zeta function. In particular, d(s, 0) has the same singularities as (s 
in {9fJe s > 0}. 
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We focus our attention in the region {9fte s > 0}. In that region, the function %^(s, 0) has 
the same singularities as 

-(s-l)-L-, 



This shows that the singularities of Z 2 appear when L s $ has 1 as an eigenvalue. In the 
next paragraph, we show that this occurs for s = | ± ir n . Then we identify the residues. 

Remark 7.4. In our conventions, r n > and we have defined the boundary values Ti Tn using 
this choice of convention. For simplicity, we will restrict our attention to s = s n = \ + ir n , 
but the analysis at s — (1 — s n ) — | — ir n would be similar. 

7.3.1. Location of the poles. For s = s n , one can check directly that 1 is in the spectrum of 
L S fi\ the eigenspace is spanned by the functionals 

/ 71n \ ir n / \ dTi r (Zi) 

(7-10) xn){zuz 2 ) = - — , 

where T irn are the boundary values of eigenfunctions of the laplacian. 

The verification is easy : let ip{ z i) be an arbitrary analytic function. What we need to 
check is that 

dT irn (zi) ^ , /2 f i v /2 dT irn (o;i) 



-/ *1 _ ^2 .TTL , 



fc_i,fc 

The last expression can be transformed to 

E 9' k M) 1+Sn [ , VW^(fe*i) w Y dT ^zi) 

k-x,k J 9 k ^K k _ 1 |£ 



— g ko z2\ 



z\ - g ko z 2 \ 



E 9k (z 2 y^ [ ^g^k^y^M)^ 



Remark 7.5. Similarly, the functionals 

dT ir (z 2 ) 



ipizt)-, 7-dzi- 

\Zi - Z 2 \ Sn 



(7.11) X l ^(zi,z 2 ) 

are eigenvectors for the adjoint L* n . 

Conversely, we need to know that 1 is in the spectrum of L s $ only if s is one of the s n ; 
and that the multiplicity of 1 is exactly the multiplicity of s n (l — s n ) in the spectrum of the 
laplacian. To prove this, we can use the relation with the Selberg zeta-function (Remark 
17. 3|) . For the latter we know that the zeros occur when s(l — s) is in the spectrum of the 
laplacian, with the same multiplicity. 

It follows also from the formula ()7.10|) relating eigenfunctions of the laplacian to eigen- 
functions of L Sn , that the characteristic space for the eigenvalue 1 of L Sn ^ corresponds to 
a diagonal matrix (there are no Jordan blocks). 
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7.3.2. The residues. We are interested in the singularities of Zi in {Jte s > 0}, or equivalently 
in the singularities of 

d z d(s,0) ^d z det{l-B (i) ' 



s,0) 



rf(s,0) Y de*(l-S a)0; 

From the previous paragraph, we know that d(s, 0) = if and only if s — | ± ir n (1/4 + r\ 
being an eigenvalue of the laplacian). For some i, the operator Bg has 1 as an eigenvalue, 
and its multiplicity rrii is the same as the multiplicity of 1/4 + r\ in the spectrum of the 
laplacian. As in the previous paragraph, we treat the case of s = s n = | + ir n ; the case of 
s = | — ir n would be similar except for the choice of a different convention in the definition 
of boundary values. We will see that the singularity of the function (|7.12|) at s = s n is a 
pole; the residue must then be given by 

d z d^- l det{l-B%) 

rrii 



8pdet(l-B®) 



Theorem 11.31 will then follow directly from: 
Proposition 7.2. 

d z d™^det (1 - b£ ) _ / X a dPS v 



s "> : —— " : Yl 



d?det(l-B?l ) r j^ rn !xdPS irj ■ 

Proof. If 1/4 + r\ is an eigenvalue of the laplacian of multiplicity m, we know that 1 is an 
eigenvalue of L Sn ^ of multiplicity m. We also know - and this is rather important - that the 
eigenvalue 1 corresponds to a diagonal block for L Sn>0 . 

Let V C C be a neighbourhood of 1 that does not meet the rest of the spectrum of L Snt0 . 
Let P Sj2 be the spectral projector on V for the operator L S 2 . As before, denote L Sj2 P Sj2 = B SjZ . 
Because we have a diagonal block, B Snt0 = P Sw ,o- Using the previous notations, for some i, 
the operator B Sn is one of the -B^o s ; it has 1 as an eigenvalue with multiplicity rrii — m - 

In the tensor product H Am , the projector P^ is of rank 1. Let V Sni o G H Am be the 
associated eigenvector; it also belongs to Ker(I — L Snt0 ) Am . By perturbation theory, we can 
find a family V sz , depending analytically on (s,z) in a neighbourhood of (s n ,0), such that 

We have 

(7.13) (I - L s , z ) Am V s , z = Xs, z V s , z 

with X s z = det(I — B sz ). 

Similarly there is a family T s z in the dual 7i* Am , depending analytically on the parameters, 
such that 

(7.14) (/ - L; tg y m z, z = \ s , z z, z 

Differentiating I7TT31 once with respect to the parameters, and applying % nj o, we get 
(7.15) 

(dX Snt o)(V Snt0 , r Snt0 )+X Sni o(dV Snt0 , T Sn ,o) = ((d(I-L Sn>0 ) Am )V Sn , , T Snt0 ) + ((I-L Snt0 ) Am )dV Sn , , T Sn , ) 
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Because (/ — L* n ) Am 7s n ,o = ^s n ,o% n ,o, the second term on each side of 17. 151 are equal, and 
17.151 amounts to 

(7-16) (dA Sn ,o)(V Sn ,o,T Sni0 > = ((d(I - L s „, ) Am )V s „,o,r Sni0 ) 

This last term vanishes if m > 1, and thus we see that d\ Snt o = 0. Iterating this procedure, 
we see that any derivative of order < m of X SjZ vanishes at (s n , 0). 

This proves, in particular, that the singularity of the function 17. 121 at s n is at most a pole, 
and that the residue we are interested in is 

as announced earlier. 

The same procedure already applied (differentiate I7.1H| then apply Zs„,o) gives after m 
steps: 



(7.i7) (a^r 1 A Sn ,o)(v Sni o,r Sni0 ) = ((^ar 1 (/-^„ ) o) Am )v s „,o,r Sni0 ) 

m— 1 

= (_l)»( m - 1)! ^((<9 s L) Afe A d z L A (<9 s L) Am ~ 1 ~ fe V Sni o, r Snj0 ) 

k=0 

The terms where L is not been differentiated disappear, because 1 — L* n vanishes on T Sn fl. 
Similarly, 

9TK,o = (-l) m m!((9 s L) Am -V s „,o,T s „, ) 

We note that d z L = L o M4 (where Ma denotes multiplication by A) and d s L = L o M T . 
Remembering that L* n T Sn fi = % n ,o we can rewrite the last two expressions as 

m— 1 

d^Kfi = {-l) m {m - 1)! ^<M T Afe A Ma A M^-^V^o, T s „, > 

fc=0 

and 

STK.fi = (-l) m m!(Mr- 1 V s „,o,T s „, ). 
Now, we can choose to write T Sn fi as 

T A ir J 

^s„,0 — /Vj=r„X(i) 

and V Snj0 as 



where X(j)>X(r) are associated to T irj by the formulae (|7.1Uj) . (|7.11|) . 

For Tj = r k = r„, we have 

^'.«S)°/'M T ' i l ?y 1 



{K1){Z 2 ) = /Xo(Sn) (0irj, 0ir fc ), 
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by the formulae of Part 5 (which could as well be applied for two different eigenfunctions 
of the same eigenvalue). Because the basis (^» r .) is orthonormal, this coefficient vanishes 
except for j = k. 
Similarly, 

and if j — k this is exactly the Patterson- Sullivan distribution applied to a. 
We finally find the expression of the residue. 

r/rn w ,T ir (dz l )T ir Adz2) 

Tfl = 7 

« s "' u rj=r n J {/•Cl){Zi, Z 2 ) | 2l _ 22 p s „ 

which is what we expected in terms of Patteron-Sullivan distributions. 

□ 

8. Classical Selberg trace formalism 

We now begin the Selberg trace formalism proof of Theorem 11.31 To prepare for the proof, 
we review the standard theory of the Selberg zeta function and trace formula and then give 
a non-standard proof which will be generalized in the next section. 

As above, we denote by {<fiir k } an orthonormal basis of A-eigenfunctions on T\G/K, with 
associated eigenvalues = 5^(1 — Sfc) with S/» = \ + irk- In particular the trivial eigenvalue 
A = corresponds to s = 0, 1 and r = ±|. 

8.1. Standard Selberg zeta function. We now review the analytic continuation and polar 
analysis of the Selberg zeta function. We refer to [V] for background. 
The Selbert zeta function is defined by 

z( s ) = n {P} n^ =0 (i-iv(P)- s - fe ), Mes>i 

where {P} runs over conjugacy classes of primitive hyperbolic elements and where N(P) = 
e Lp where Lp is the length of the corresponding geodesic. 

The logarithmic derivative of the Selberg zeta function s _^ 2 ^( s ) is defined for 3fte s > 1 
by the formula (see jVj, (5.1.5)) 

v ' ; s- 1/2 Z K ' ^sinhL 7 /2 

In this formula, we sum over all closed orbits 7 of the geodesic flow and L 1 is the (positive) 
length of 7. 

Z'(s) 
Z(s) 

together with the 'trivial poles' at s = —k, k — 0, 1, 2, 3, . . . 



Theorem 8.1. admits a meromorphic continuation to C with poles at the points s = s n 



Proof. We review a few features of the standard proof to draw attention to some impor- 
tant technical issues which might be confusing for the more general versions to come. By 
definition, we have for 9?e s > 1, 

1 Z'(s) _ 1 Z'(a) _ L l0 r 1 -( S -1/2)L 7 L_p-(a-l/2)i 7 \ 

s-1/2 Z{s) a-1/2 Z(a) ^7 sinh L 7 /2 L 2s-l 2a-l J" 
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See for instance [EJ, (5.1.5). 

To analytically continue the formula, one applies the Selberg trace formula (cf. jVj Theo- 
rem 5.5.1) with the test function 

i / 1 2 \ 1 1 

n{- + r ; s; a) 



•4 ' ' ' ( s _i)2 + r 2 ( a _l)2 +r 2- 

The Fourier transform of h(j + r 2 ; s; a) is 

^ s; a ) = _J_ e -(^i)M _ _J_ e -(a-§)M_ 

yv ' ' ; 2s -1 2a -1 

We note that the rate of decay of h{\ + r 2 ; s; a) as r — > oo reflects the singularity of \u\ at 
u — 0. In the case of a smooth compact quotient, the result is (see jV] Theorem 5.1.1; see 
also |Salj ) 

1 Z'(s) l Z» Vol(r\G) y^oo / 1 1 \ 

Z(a) (a-i) Z(a) tt 2^fc=0 U+fc a+*J 

(8.1) 

I V"^00 / 1 1 

T 2^n=0 ^ (s _i)2 +r 2 (a _i)2 +r 2 

We note that the eigenvalue series on the right side would diverge if we only used the 
formula for ^j^j, but it converges (away from poles) if we subtract or take one derivative. 

These formulae give a meromorphic continuation of to C and show that the poles 
occur at values of s for which there exists an eigenvalue A n satisfying A n = s(l — s), or at 
negative integers. 

8.2. Convolution operator approach. As sketched above, the Selberg trace formula in- 
volves a Fourier transform duality. We will need a more group theoretic approach for the 
generalizations in the next section, namely the approach in |GGPj to the Selberg trace for- 
mula as a formula for the trace of the convolution operator corresponding to a i^-bi-invariant 
function \- 

We denote by Sq } q(G) the continuous functions satisfying x(kigk 2 ) = 4>(g) for all k\, k 2 G 
K. The associated convolution operator is defined by 



R x= x(g)R g dg, 

JG 



where R g f(x) = f(xg). There exists a unique (up to scalars) eigenfunction \l/ s of Q of 
eigenvalue s(l — s) in S^q. The spherical transform S : C£°(G) fl S^q — > PW m is defined by 



Sf(s) = [ f(g)$ s (g)dg. 

JG 



' G 

Its range is the subspace of the Paley- Wiener space 

(8.2) PW(C) = {f e 0(C) :3k\/N>0: \f(x + iy)\< Ce fcN (l + \y\)~ N } 

with a certain symmetry which we will not need to recall here (see p. 31). Here, 0(C) 
denotes the holomorphic functions on C. 

We also denote the Mellin transform M : C^(A) -> PW(C) by 

Mf(s)= r f (a)a s -, 
.In a> 
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where we identify f(a) as a function of the top diagonal entry of a. Note the non-standard 
sign of the exponent, which is chosen to be consistent with jLj |Z]. 

The basic Selberg trace formula for a smooth compact quotient (in the form stated in 
[GGP]) states that 

(8.3) V S X (2ir k ) = Vol{V\G) X {e) + V / xGrSs)^, X e S , 

7~, Jg^\g 



r k {7} 



where the sum runs over the principal and complementary series representations (counted 
with multiplicity), where G 7 is the centralizer of 7 in G (similarly for V) For 7 7^ e, T 1 \G 1 
is a closed geodesic. 

The orbital integral on the right side of (|8.3|) may be expressed in terms of the so-called 
Harish-Chandra transform as follows: If X $00, there exists x° on D = G/K such that 
x{d) = X D {9 ' 0) where x D (re i61 ) = x D { r )- I n the proof of [Zj, Proposition 2.6, it is shown 
that x( n u lan u) — X D (\^ft;\)i wrfc h uj = ■ With some routine manipulation ([Z], pages 
55-56), we get 



(M H X (a) = \ a -a^\\uj\ j^ X D {{l + ^i) 1,2 )d 



u. 



and thus 

(8.5) / X{g- l 19)d9 = V ° l{V ^ H X {a^ 

JG^\G l a 7 ~~ a 7 I 

where ^oZ(r 7 \G 7 ) is the length of the closed geodesic. We further have Sx = MHx, so we 
finally obtain 

(8.6) MH x(2*r k ) = Vol(T\G)x D (0) + £ — # xK)- 

»•* {7} 7/ 

This approach leads most naturally to the zeta function 

(8.7) K(s;1): ^(_^_) (c08hV2) -,.-V, 

In generalizations to non-constant automorphic forms, we begin with (|8.7j) and then relate it 
(and its generalizations ()9.5|) to non-constant automorphic forms) to the usual zeta functions 
Z(s: a). 



9. Dynamical zeta functions and Selberg trace formalism 



This section is concerned with the zeta functions 




(Ke s > 1). 



Theorem 9.1. For each automorphic form o = <p irk , X + <f) irk ,tfj m , Z(s;a) is absolutely con- 
vergent in 3fte s > 1 and admits a meromorphic continuation to C . Except for the trivial 
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representation a = 1, the only poles in 3?e s > occur at values s = | + ir for which | + r : 
zs an eigenvalue of A, and the residue is given by 

Res s= i +ir Z(s; a) = /i ( \ + ir) Y] [ adPS irj . 

>r\G 



j: r ?= 



This proves a special case of Theorem 11.31 in which the function a has components in 
a finite number of irreducible representations. We briefly sketch the extension to analytic 
symbols in the final section. 

The proofs are based on a generalized Selberg trace formula introduced in [Z] for the 
traces TraR x on L 2 (T\G) of the composition of R x with multiplication by a. Here, a is a 
Casimir eigenfunction of weight m and R x is a convolution operator with kernel X £ £>m,n(G), 
where 5 TO>n denotes the functions xid) on G satisfying x(k^gkg) == e tm ^ e m6 'x(g) , where 
kg = expOW G if. The eigenspaces of fl on S m ^ n {G) are one-dimensional, spanned by the 
spherical function $ m . njS of ^-eigenvalue s(l — s). We will only be considering the case n = 0, 
and denote the associated normalized spherical function by $ m>s . Our normalization follows 
[HI IZ]- The spherical transform S m : C£°(G) fl S m fi — > PW m is defined by 

S m f{s) = / f{g)$- m ,s{9)dg- 
Jg 

Its range is the subspace of the Paley- Wiener space (|8.2|) with a symmetry depending on m 
which we will not need to recall here (see [ZJ, p. 31). 

We will also need a variety of Harish- Chandra transforms which depend on the weight m 
and also on the type of representation P ir .,X>+. There is a canonical one, defined as follows: 
Let x £ 5m,o and let (see [ZJ page 57 for (i) and page 49 for (ii)): 

(9-2) H mX (a) = la-a-^JZx^m^ 0" = S£) 

Here, if x G S^o then there exists x D on D = such that x(g) = x D (fl' ' 0) where 

X D (re ie ) = e 1 ^ X n(r). In the proof of [Z], Proposition 2.6, it is shown that 

D( u+i \ _ j!Z$(a,u) D (\ u+i \ \ 



x(n- 1 an u ) = x D (^t) = e^^ X D (\i 



-IUJ I 



with 

(9.3) e«"> = <ft±^M±^l , where , = £±£l, „ = 

(it + iuj)(u 2 + l) 1 /^ a — a 1 

With some routine manipulation (see [ZJ, pages 55-56), we get 

(9.4) H mX (v) = \a-a-'\\u\v m ^ J R (l + ^)- m V((l + ^) 1/2 )d 



u. 



From the Selberg trace formalism viewpoint, it turns out to be most natural to work first 
with auxiliary dynamical zeta-functions 7Z(s; a) that do not seem to arise in the thermody- 
namic formalism. When a m has weight m we put 

(9.5) K{s- a m ) := V I ) (tanhL 7 /2) m/2 (coshL^)- 2 ^ 1 / 2 ). 

\ smh L 7/2 / 
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We then express Z(s; a) in terms of TZ(s; er) to obtain results on the analytic continuation 
of the latter. This somewhat circuitous route comes about because the trace formula is on 
the 'quantum level' and therefore does not quite produce the 'classical' zeta-function. 



9.1. Forms of weight in Vi r . In this section, we prove Theorem 11.31 for the case a = 
Z(s;<j>ir k ). 

In this case the auxiliary zeta function has the form 
(9.6) K{s;M := £ (J^^J (coshL 7 /2)- 2 ^). 

Theorem 9.2. lZ(s; <f>i rk ) admits a meromorphic continuation to C with poles at s = \ + 
%Tj — k, k = 0, 1, 2, . . . , where | + r| is an eigenvalue of A, and with 



Res« 



<t>ir k ) = M* r + n) \^ir k ,PS rj ) ■ 



j: r 2 .=r 



Proof. We assume throughout that 4>i Tk J_ 1, so that the identity term on the 5^ 7gr side 
of the trace formula vanishes and so that the trivial representation term with r — \ also 
vanishes. After the proof, we remark on the case 0j rfe = 1. 

By Proposition 2.12 of [Z] (applied in the continuous series case), we have 

(9.7) BPK^, <p irn )MH oX (2tr n ) = £ -^jL_ i^ fcX (a 7 ), 

where Hq is defined by (I9.2|) - ()9.4|) . and where (see [Z] page 57 for (i) and page 49 for (ii)): 

(9-8) m rkX (a) := \a - a^\ JZ o F irk>0 (^) X D \^t)du. 

Here, a = e L l 2 and ^ir fc , o(zif) * s defined in (|5.7jl . We note that the identity term on the right 
side vanishes by orthogonality. 

Remark 9.1. (i) We note that we do not use Proposition 2.10 of [Z], which gives a less 
convenient zeta function. Although Proposition 2.12 of [Z. is only stated for symbols in 
the discrete series, it is valid for the continuous series as long as we use the corresponding 
expressions (given in [Zj Corollary 2.4) for the integrals I 7 (er)(n u ) in jZ] (2.2) . 

(ii) A priori, the right side of ()9.8j) should also contain the term 

\a-a \ G irkfl (-—-)x ( — -—)du, 

J-oo -2i U + IUJ 

but x° is a radial function since it has weight zero and X D (|^p^;|) is even in u while 
Ginfiiz^) is odd. Hence this integral vanishes (cf. Proposition 2.7 of jZ]). 

By (jOjl . we have 



(9.9) 



H oX (v) = \a-a-*\\u;\ f RX D ((l + ^) l/ y U , 
Hi irh x(v) = {a-a-^ujl J R F irkfi (^ X D ((l + ^) 1/2 )du. 
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We now define Xs(g) £ 5*0,0 by the rule that 

Xf(r) := V - - f { , 0<r<l. 
Using ()9.9)1 and the fact that \a — a _1 ||u;| = (a + a" 1 ), we obtain 

TC 



(9.10) 



(ii)H QXs (a)= = Jgfc (a + a-r 2(s - 1/2) 



If we substitute \s into the right side of the trace formula ()9.7|1 . we obtain the desired zeta- 
function TZ(s; <pir k )- Therefore, the left side of the trace formula ()9.19j) gives a meromorphic 
continuation of 7Z(s; <pir k )- By Theorem 11.31 we have 

(Op((j) irk )(j) ir ,(j) ir } = (<j) irk ,PS ir }{i c irk (l + ir), 

hence 

(9.11) n{s-4> irk ) = J2(&r k ,PSir n )ti rk (~+ir n ) MH oXs (2tr n ). 

n=0 

By flmHii)) we have 

MH oXs (2zr) = ^Lf^a^ia + a-iy^/V^ 

= [°° e 2iri (cosht)- 2(s - 1/2) dt 

Ms} r( S -(|+ir))r(s-(|-ir)) 
/*?r fc W T(2 S -1) 

For the last line we refer to |Zj (p. 60). 

In conclusion, we obtain (at least formally) 

(9.12) K(si M = f>, rt , p S „., /oW^(^^r (a -(^^))r(.-(^-^)) 

„=0 PirfeWJ 1 l ZS -U 

We note that — ^ 2 — = B(s — (| + ir), s — (| — ir)). As above, we assume that 

0i rfc _L 1, so that the trivial representation term vanishes. Regarding the convergence of the 
right side, we note that by (|5.fij) and ()2.10j) . as \r n \ — > oo, 

c (\ i • \ -i/2 

r(a - (| + ir n ))r(a - (§ - ir n )) ~ e~f (l^+^l+l^-^l)| r „ + Ssl"*"- 1 ^ - Ofsl"^- 1 . 

Since ((p irk , PS irn ) = O rk (l) as n — > oo, it follows that the series converges absolutely in the 
critical strip away from the poles and defines a meromorphic function. 
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There are simple poles at s — \ ± ir n where | + r 2 is an eigenvalue of A. In the case 
where the multiplicity of the eigenvalue equals one, the residue at s — \ + ir n equals 

{ct> irk ,ps irn f° 2 llj{lV n )) ^ rSlr^ = ^ + ^)<</w^ r j, 

as stated. In the case of a multiple eigenvalue one sums over an orthonormal basis of the 
eigenspace. 

□ 

9.2. Z(s; 4>ir k )- Now we deduce properties of Z(s; 4>ir k ) from those of 1Z(s; 4>ir k )- 
We introduce the measure 

We note that 

Z{s- ct> irk ) = J °° e-^ L dQ{L-<t> irk ), 

-(•-3) 



(9.13) 



Lemma 9.3. We /iave: 

oo 

Z(s; ir J = ^ £ m (s, n)7£(s + n; <f> irie ), 



n=0 



™ i.. t. — n-t. j I. — « \ / \ ± / 



1 \ / I 

2 



m,fci,...,fc m =0;fci+---fc m =n 

Proof. By elementary manipulation, we have 

(9.14) Z{s^ irk )= / (l + e-^^de^js;^), 

./o 

where 

de (^ 8 ;^) = E( s ^^} 2 ) (( cosh V2)T (s4) 5(l-l 7 ). 

We then change variables to y = (coshL/2) 2 , and note that e~ L / 2 = ^/y — -y/y — 1 to obtain, 

( 9 - 15 ) , , 



where 
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By repeated use of the binomial theorem, there exist coefficients B m (s,n) such that 

25-1 oo 



(l-Jl^) =y-^Y.B m {s,n)y 



Cancelling the factors of 2/ ±( - 2s x \ we thus have 

Z{s,4>ir k ) = J2n=o B m(s,n) J^y- n d^(y;s;(f) irk ) 

(9.16) 

Since the poles of 1Z(s + n, 4>ir k ) are the shifts by — n of the poles of 1Z(s, 4>ir k )i an d since 
the non-trivial poles of TZ(s, ip m ) in 3?e s > lie only at the points s = | + ir, only the term 
n = in the series contributes non-trivial poles to the critical strip. 

Writing out B m (s, 0) as a sum ^ mfcl u =0 . k + ... k =n , we see that the only term has 
m = = kj (for all j = 1, . . . , m). Thus, 



#es s= i +ir Z(s; 0irj = Res s= i +ir K(s] Yll , 
This completes the proof of Theorem 19.11 in the case a = <pi rk 



□ 
□ 



Remark 9.2. As a check on the calculations, we observe that in the case <fri rk = 1, (4>ir k , PSi Tn ) - 
1 for all n, yUj rfc (s) = fio( s ) an d we get 

(9.17) 

+ 2^n=0 r(2s-l) 

The series converges rapidly to a meromorphic function with simple poles at s = | ± r n — k 

(k = 0, 1, 2, ... ), the residue at s = \ ± r n — k being . Thus, Lemma 19.31 shows that 
Z(s; 1) has simple poles in the critical strip with residues equal to 1. The formula (|9.17|) also 
follows from the standard (Fourier transform duality) Selberg trace formula Theorem 
4.3.6) by using the integral formula (cf. [WW], Exercise 24) 

. 1 1 . 1 f cos(2ir n u)du 

Bis — - — zr n , s — - + ir Tl 



o o 1 " nj — i ; i2s-l/ \ 

^ 2 4 2 Jr cosh (u) 

and the fact noted above that — ^^(^i) ~ ~ = B(s — (| + ir), s — (| — ir)). 

9.3. Forms of weight ±2 in TV. In this case both sides of the trace formula equal zero 
due to time reversibility. By Propositions 12.51 and 13. 3| each side of the trace formula equals 

/[ X+6 ir ,ds\ {L-l X +<f>ir k ds \ 

zero, noting that ( ^J^ ) + = 0. 
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9.4. Weight m in X>^. We now prove Theorem 19 .ll for a = ip m G Z>+. The anti-holomorphic 
discrete series case is simply the complex conjugate and is omitted. 

The proof is similar to the case 2(s; <pir k ) but involves the higher weight analogue zeta- 
function: 



(9.18) n(s-^ m ) :=y( ^™ dS \ (tanhL 7 /2) m/2 (cosh^^)-^- 1 / 2 ). 

We begin the proof with an analysis of its meromorphic continuation. 
9.4.1. Meromorphic continuation of lZ(s; ip m ) . 

Theorem 9.4. TZ(s;ip m ) admits a meromorphic continuation to C. In the critical strip, 
its poles occur at s = | + ir such that j + r 2 is an eigenvalue of A , with residue /io(| + 



Proof. We study lZ(s;ip m ) using the trace formula given in j7j, Proposition 2.12: 

sinh L 7 / 2 



(9.19) £)<0p(«^, <f>ir n )MH mX (2ir n ) = ]T ( .^"^ ) <*K,- 



n=0 7 

where (see jZ] page 57 for (i) and page 49 for (ii)): 

(i) H mX (a) = la-a-'lj^X^^du, 

(ii) Hi X {a) = \a-a^\^ oo {u + i y m /\ D {^t)du. 



(9.20) 



We caution that in the definition of (|9.2Uj) (ii) we follow a slightly different notation 
convention in [7] whereby we multiply the integral by \a — a~ 1 \ as for H m 

The integral uses the notation of (|9.2|) - (j9.3|) . We simplify the expressions in (|9.2(J|) by 
further using these identities to obtain (see also [Zj, pages 55-56) 
(9.21) 

(i) H mX {v) ~ --HI..I.W4 f M , v—l \~m/A^,D ( (-[ , „-lW 2 > 



a 


-a" 1 ! 


a; 


a 




a; 



(ii) = la-a- 1 !^!^/ 4 J R (« + i)- m/2 (l + ^)- m V((l + ^) 1/2 )^ 

We now define Xs{g) G £m,o by the rule that 

r- m/2 xf (r) := 1 d , J , < r < 1, 
where (see jZ], Proposition 3.6) 

(-i) m / 2 7r2 2s+2 - m / 2 r(-2s + 2) 



Mr, 



( S ) = / ( M + 2)- m / 2 ( M 2 + l)-^M 



r -(2s + i-2)r(- s )r(- s + ffi)- 
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_! W2 



Since (1 — v) s = (a + a" 1 ) -25 ', and \a — a _1 | |co>|t> m / 4 = (a + a -1 ) ( j , we have 

y / \ m/2 

(<) ^X.(o) = (t^) (a + a^)- 2 ^), 



(9.22) 



It follows first that if we substitute \s into the right side of the trace formula (|9.19|) is the 
desired zeta-function 1Z(s;ip m ). Therefore, the left side of the trace formula f)9. 19j) gives a 
meromorphic continuation of 7Z(s; ip m ). By Theorem 11.31 we have 

hence 

(9.23) K(s; = ^(W, PS irn )i2 d m (- + %r n ) MH mXs (2ir n ). 

n=0 

By dS^ ii)) we have 

MH mXs (2ir) = ^r« 2ir (^) m/2 (« + 0- 2(s - 1/2) ? 

Putting things together, we obtain the discrete series analogue of (j9.12|) . 
(9.24) 

1 / \ /*oo 

K(s; rj, m ) = V> m , PS irn )» d m (- + ir n ) J^l / e 2 ^ (tanht) m/2 (cosht)- 2 ^ 1 / 2 ^. 

The integral is more complicated than its zero weight analogue, but as tanht = 1 + r(t) with 
r(t) = 0(e~ 2 '*!), we may write 

/oo poo 
e 2irnt (tanht) m/2 (coshty^-^dt = / e 2irnt (cosht)- 2is - 1/2) dt + R 2 (s, r n ), 
-oo J — oo 

where 

/oo 
e 2^t r ^( cosht )-2( S -l/2) df _ 
-oo 

The first term of ()9.25j) gives the expression in the weight zero case analyzed above. Hence, 
the sum over r n with this term converges, and the poles and residues of TZ(s; ip m ) on Se s = 
1/2 due to this term are the same as for 

Ms)»t(h+ir n ) r(s-(i+ir n )r(s-(|-tr„)) 



&(») r(2s-l) 

Ho{\-\-ir n ))^{\-\-ir n ) 



There are simple poles at s = ~ + ir n and the residue is (ip m , PSi r )- — — : — ■ — 

(■0m 3 PSir n )^o{\ + iy n ). Summing over an orthonormal basis of lowest weight vectors of 
gives the stated expression. 
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To complete the proof, it is only necessary to observe that the second integral i?2(s, r n ) is 
holomorphic in the region 3?s > — ~. It is also rapidly decaying in r n . Therefore it does not 
contribute any poles or residues to TZ(s; ip m ) in the critical strip. 



9.5. Z(s; i/) m ). Now we deduce properties of Z(s; ipm) from those of 7Z(s; ip r , 
As in the weight zero case, we introduce the measure 

dQ(L^ m ) = {^^J2 ) ( tanh V2) m/2 6(L-L 7 ). 



We note that 
(9.27) 



Z(s-M = / °° e ' is ^ )L (tanhL/2)- m/2 dQ(L;^ r , 

-(*-§) 



Because the factor (tanhL/2)~ m//2 is somewhat inconvenient, we also consider 

= (tanhL 7 / 2 y 



(9.28) 

Lemma 9.5. We have: 



II e -( s -I)L 7 _ 



Z(s; ip m ) = ^ B m(s, n)K{s + n; ip m ), 



n=0 



where B m (s,n) is the same as in Lemma \9. 'J[ 

Proof. We use similar manipulations as in the weight zero case. We now have 

roc 

(9.29) Z(s^ m )= / (l + e- L ) 2s ~ l dQ(L-s-il> m ), 

Jo 

where 

dQ(L,s;ij m ) = J2 (| s ^ / L m /2 ) ( tanh V2) m/2 ((coshL 7 /2) 2 )- (s -^ 5(L - L, 

We change variables as before to y = (coshL/2) 2 , and obtain as in (J9.15)) . 

(9.30) Z{sM = ^{2yf^( K l-^l^l) 2S ' l d^{y-s^ m ) , 
where 

d*(y;s-M=Y,( ^f/o ) (tanhL 7 /2) m / 2 (y,)^ 8{y-y,). 



7 



sinh L 7 /2 

As in the weight zero case, we then have 

Z(s,ip m ) = Y.n=o B rn(s,n) y- n d^(y]S]i). 

(9.31) 



□ 
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Since the poles of 1Z(s + n,ip m ) are the shifts by — n of the poles of 7Z(s,ip m ), and since 
the non-trivial poles of 7Z(s, ip m ) in 3?e s > lie only at the points s = \ + ir, only the term 
n = in the series contributes non-trivial poles to the critical strip, and as above this term 
has m = = kj (for all j = 1, . . . , m). Thus, 



Res, 



j+ir 



Res c 



j+ir 



□ 



To complete the proof of the theorem, we now observe that 

L ^rnds 



(9.32 



Since 



(9.33) 




sinh L 7 /2 
and since 

(_ s _l+i)L 



1 - (tanh LJ2) 



m/2 



< oo, 3?e s > 0, 



by the prime geodesic theorem, it follows that Z(s,ip m ) has the same poles and residues in 
the critical strip as Z(ip m ). 

This completes the proof of Theorem | 



9.6. The analytic case. We briefly dicuss the extension of Theorem 19.11 to more general 
analytic symbols o G C w (r\G). 

By Proposition 12.31 by a similar calculation as in Corollary 16.161 we have 



(9.34) 



Z(s-a) 



E 



— -Z(s; (p ir . 



(a 



(3Js > 1) 



Here, we interchanged the summation over 7 and over which is justified by Proposition l2.4l 
and the prime geodesic theorem. The significant new question is the analytic continuation 
of (J9.34|) to the critical strip. If a has only finitely many components, then of course the 
analytic continuation of the sums comes down to that of the individual terms. To extend 
the analysis to infinite series, we recall that the coefficients 

(9 35) — ^ ir ^ ^m) 

are rapidly decaying in rj, resp. m (see (|6.7J) ) as long as a G C°°. Although we will not 
investigate it here, for a G we expect that the coefficients are of exponential decay. 

To extend the analysis to infinite series, it is necessary to have estimates of the zeta 
functions Z(s; 4>i r .), resp. Z(s; ip m ) in the parameters Tj, resp. m. By Lemma l9~31 and ()9.31|) . 
the estimates comes down to the estimates of 7Z(s; (fi ir .),7l(s] ip m ) in the r k , m parameters in 
the critical strip. By the trace formulae, this in turn reduces to estimates of the coefficients 



(9.36) 



ir k \ 2 



+ ir r . 



II' 



in the joint (r k ,r n ) resp. (r fc ,m) parameters. Proposition 15.31 reduces these estimates to 
ratios of T functions. Hence, one can obtain the analytic continuation properties of Z(s; a) 
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as long as the coefficients (|9.35j) decay rapidly enough to offset the growth of the T factors 
in ([9.36)1 . But we will not determine the optimal statement here. 
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